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Abstract 


An  extended  Kalman  filter  is  used  to  predict  a  kinetic  kill  zone  for  use  in 
aircraft  self-defense  versus  homing  missiles.  The  analysis  is  limited  to  an  in-the- 
plane  analysis  and  focuses  on  finding  the  model  parameters  which  have  the 
largest  impact  on  the  predicted  kill  zone.  No  attempt  is  made  to  optimize  the  design 
of  the  filter  model  itself.  The  analysis  computes  the  kill  zone  relatitve  to  an 
assumed  aircraft  trajectory  using  strictly  filter  computed  statistics.  No  Monte-Carlo 
simulations  are  used  throughout  the  thesis.  The  filter  assumed  to  be  on  the  evading 
aircraft,  uses  an  onboard  laser  radar  (ladar)  to  provide  measurements  of  aircraft- 
to-missile  relative  range,  range-rate,  line-of-sight  and  line-of-sight  rate.  The 
missile  is  assumed  to  be  in  a  post  burnout  coast-to-intercept  phase  of  flight. 
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FEASIBILITY  ANALYSIS  FOR  PREDICTING  A  KINETIC  Kill  ZONE 
FOR  AIRCRAFT  HOMING  MISSILE  DEFENSE 


I  Introduction 

1.1  Background 

The  class  of  munitions  called  "expendables",  chaff  and  flares  has  always 
been  a  primary  source  of  endgame  aircraft  self  defense,  primarily  because  they 
are  economical.  The  economic  advantage  comes  from  being  able  to  defend  against 
a  large  class  of  missile  types  with  only  a  few  standard  expendables. 

Due  to  increasing  complexities  in  missile  seeker  Electronic  Counter- 
Countermeasure  (ECCM)  technologies,  seif  defense  expendables  used  for  Electronic 
Countermeasures  (ECM)  purposes  have  become  increasingly  more  diverse  and 
specialized,  limiting  the  expendables'  generic  usefulness.  Optical  and  electro- 
optical  missile  seekers  have  become  common,  and  their  sophistication  is 
demanding  new  expendable  technology  in  the  visual  and  Infrared  (IR)  spectrums. 
For  example,  IR  missile  seekers  are  using  spectral  color  discrimination  and  source 
rise  times  to  help  reject  IR  flares.  Radar  based  seekers  employ  home-on-jam, 
special  frequency  filtering  and  pulse  coding  techniques  to  defeat  radar  ECM 
techniques.  The  spread  of  ECCM  technologies  forces  ECM  technologies  to  diversify 
to  provide  adequate  multispectral  coverage,  thus  reducing  the  economy  of  using 
expendable  munitions  for  aircraft  self  defense. 


An  expendable  mini-ball  shot  cannister  might  potentially  restore  the 
generic  usefulness  inherent  in  expendable  munitions  [20|.  Conceptually,  such  a 
kinetic  kill  cannister  would  allow  one  expendable  package  to  be  deployed  in 
virtually  any  endgame  missile  scenario  because  the  kinetic  kill  cannister  does  not 
defeat  the  missile  by  jamming  the  missile's  seeker  head,  but  rather  by  hitting  the 
missile  as  depicted  in  Figure  1-1.  The  successful  employment  of  a  kinetic  energy 
impact  device  would  negate  current  seeker  ECCM  capabilities  provided  by  the 
proliferation  of  missile  seeker  technologies.  However,  in  order  to  use  such  a 
defense  mechanism,  one  must  accurately  predict  the  future  flight  path  of  the 
hostile  tracking  missile.  This  prediction  is  made  using  noise-corrupted  and 
incomplete  measurements  from  the  evading  aircraft's  sensors.  The  predicted  final 
time  of  the  kinetic  intercept,  must  exceed  the  sum  of  the  kill  zone  estimation 
processing  time  and  the  kill  cannister  deployment  time,  and  must  still  allow 
enough  remaining  missile  flight  time  to  prevent  the  disabled  missile  from 
entering  into  a  lethal  fusing  range. 

1.2  Problem  Statement 

This  thesis  explores  the  use  of  an  extended  Kalman  Filter  to  estimate  the 
future  flight  path  of  a  homing  missile  tracking  a  friendly  aircraft.  The  missile's 
predicted  flight  path  is  estimated  using  noise  corrupted  measurements  available 
from  an  onboard  Laser  Radar  (ladar),  and  the  aircraft's  Inertial  Navigation 
Subsystem  (INS),  and  the  estimated  flight  path  is  used  to  define  a  kill  zone  volume. 
The  ladar  provides  range,  range-rate,  angle  and  angle-rate  measurements  of  the 
missile  relative  to  the  aircraft.  The  INS  provides  aircraft  velocity  measurements. 
Acceleration  measurements  are  also  assumed  to  be  available  from  a  g-meter  or 
other  source.  Aircraft  maneuvers  are  restricted  to  readily  predictable  constant 
speed,  constant  angle-rate  turns  (see  Section  1.4).  The  calculated  kill  zone  volume 
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1.  The  missile  at  B  is  tracking  the  Aircraft  at  A  to  an  intercept  point  h. 

2.  The  aircraft  is  taking  measurements  of  the  tracking  missile  and  has  reached 
steady-state  in  the  error  estimation. 

3.  While  at  A,  the  aircraft  estimates  that  the  missile  will  be  contained  within 
D,  when  the  aircraft  reaches  C. 

4.  The  aircraft  deploys  the  countermeasure  from  A,  to  intercept  D  while  the 
aircraft  flies  the  assumed  trajectory  to  C. 


Figure  1-1  Illustration  of  Hard  Kinetic  Kill  Zone 
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is  determined  by  propagating  the  filter  estimated  range  and  angle  covariance 
errors  forward  to  the  desired  intercept  time. 

The  boundaries  of  the  kill  zone  are  estimated  through  a  three  step  process. 
First,  during  the  missile  acquisition  phase,  missile  detection  and  filter  initial 
conditions  are  provided  to  the  filter  from  a  Missile  Warning  Receiver  (MWR). 
Second,  during  the  estimation  phase,  the  filter  uses  the  measurement  updates 
provided  by  the  aircraft  sensors  to  determine  the  errors  in  MWR-provided  states. 
Finally,  during  the  prediction  phase,  the  filter  propagates  the  states  forward  to  the 
desired  intercept  time.  Because  the  filter  is  predicting  forward  in  time,  no 
measurement  updates  are  provided  during  the  prediction  phase,  causing  the  filter 
covariances  to  diverge. 

The  kinetic  kill  zone  is  given  as  a  horizontal  plane  volume  cross-section 
throughout  this  study,  since  ail  models  and  analysis  techniques  used  in  this  thesis 
are  defined  in  the  horizontal  plane.  Actual  calculation  of  the  kill  zone  volume 
proceeds  as  follows.  The  desired  future  intercept  time  is  chosen.  The  filter  states 
are  propagated  forward  to  the  desired  intercept  time.  A  from  the  filter  states,  a  new 
aircraft  position  is  estimated  for  the  desired  intercept  time.  The  missile's  position  is 
then  calculated  from  the  filter's  range,  5R(t),  and  angle,  0(t)  estimates  defined 
relative  to  the  new  aircraft  position  at  the  desired  intercept  time,  and  the  filter's 
estimated  la  bounds,  3*e(tf)  and  0e(tf),  at  the  desired  intercept  time.  These 
parameters  determine  the  two  dimensional  intercept  cross-section  shown  in  Figure 
1-2.  The  volume  cross-section  representing  the  final  kill  zone  shown  in  Figure  1-2 
is  calculated  as 

Volume  Cross-Section  =  SH(tf)  6SR(tf)  50(tf)  (1-1) 

where  the  limits  of  integration, 9t  and  0,  are  defined  by 
fR  over  the  interval  5R(tf)  ±n3Mtf), 

0  over  the  interval  0(tf)  ±  n  0e(tf), 
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and 

9t(tf)  is  the  estimated  range  at  the  final  time, 

0(tf)  is  the  estimated  l.ine-of-Sight  (LOS)  angle  at  the  final  time, 
;He(tf)  and  0e(tf)  are  the  filter’s  estimated  lo  bounds, 
n  determines  the  number  of  standard  deviations  the  volume  encloses. 


Parts  of  the  problem  which  require  investigation  include:  ( 1 )  How  do  the 
estimation  errors  of  SR(t)  and  0(t)  propagate  as  a  function  of  the  time-to-intercept 
(Tjj)  parameter?  (2)  How  do  aircraft  maneuvers  affect  the  TTj,  central  intercept 

point,  and  associated  errors  which  define  the  volume  cross-section?  (3)  What  laser 
measurement  parameters  and  assumptions  have  the  greatest  impacts  on  the  ability 
to  predict  the  volume  cross-section?  (4)  How  do  dynamic  elements  used  (e.g., 
missile  dynamics,  ladar  errors,  Kalman  filter  models,  navigation  errors,  etc.)  affect 
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the  final  volume  cross-section  parameters/  (5)  How  do  scenario  initial  conditions 
such  as  initial  missile/ aircraft  state  values,  allowed  aircraft  maneuvers,  and 
specifics  of  a  given  scenario  problem  affect  the  estimation  results/ 

Figure  1-3  illustrates  the  relationships  between  the  various  timing  elements 
which  hypothetically  combine  to  define  the  Ty]  parameter  required  to  predict  the 
kinetic  kill  zone  of  a  tracking  missile.  In  Figure  1-3  time-to-go  (T^)  is  defined  as 
the  scenario  time  remaining  before  the  missile-to-aircraft  range  equals  the  missile 
fusing  range,  Ts  is  the  filter-state  propagation  time  between  measurement  updates, 
Tc  is  the  processing  time  for  the  estimation  filter's  last  measurement  update,  Tt  is 
the  filter's  required  estimation  processing  time  to  predict  a  volume  cross-section, 
and  is  the  cannister  (deployment  time  to  intercept  the  volume.  Note  that  the 
possibility  of  achieving  kinetic  kill  necessarily  requires  that 

tTT,  > TT1  =  Tp  +  Te  +  Td  ( 1-2) 

where 

Ttc  is  the  time  to  aircraft-missile  impact, 

Tji  is  the  time  to  kinetic  kill  of  the  missile, 

Tp  is  the  final  measurement  processing  time, 

Te  is  the  volume  estimation  time,  and 
Td  is  the  cannister  deployment  time. 

1.3  Summary  of  Current  Knowledge 

The  kinetic  kill  zone  estimation  problem  is  largely  a  system  level  integration 
of  several  smaller  problems.  Several  of  these  smaller  problems  have  aiready  been 
adequately  solved,  and  other  parts  are  concurrently  under  investigation  which 
include  use  of  a  guided  kinetic  kill  projectile,  and  state  estimation  of  a  homing 
missile  using  a  ladar. 

Predicting  the  future  state  of  a  homing  missile  first  requires  accurate 
estimation  of  the  missile’s  current  kinematic  states.  Cusumano  and  Deponte  used  an 


1-6 


Figure  L-3  TTI  Parameter  Flow  Chart 

(This  figure  illustrates  a  typical  timing  sequence  necessary  to  complete  a  missile  intercept.) 
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onboard  radar  to  track  a  homing  missile  {5,1 1).  Measurements  taken  by  the  radar 
were  processed  by  an  extended  Kalman  filter,  which  produced  sufficiently  accurate 
state  estimates  of  the  homing  missile  for  tracking  purposes.  An  in-depth  analysis 
of  the  Kalman  filter  estimation  problem  was  performed  to  minimize  estimation 
errors  of  the  missile  states.  Cusumano  and  Deponte  focused  on  state  reduction  of 
the  missile  dynamics  model.  The  authors  developed  a  5th  order  missile  dynamics 
truth  model,  and  tested  reduced  order  missile  models  (zero  to  4th  order)  for  relative 
performance  in  a  filter  model.  The  test  results  indicated  that  filter  performance 
increased  only  marginally  as  the  filter  model  complexity  was  increased  from  a  zero 
order  to  a  fourth  order  missile  model  [5:1.  The  kinematic  equations  of  motion 
between  the  aircraft  and  missile  had  been  developed  in  previous  thesis  work  by 
Lutter  [8,101.  Together,  the  thesis  work  by  Lutter,  and  that  by  Cusumano  and 
DePonte,  provide  a  good  foundation  for  target  tracker  extended  Kalman  filtering. 

Negast,  in  his  thesis,  analyzed  the  performance  of  various  Global  Positioning 
Satellite  (GPS)/INS  integration  approaches  versus  the  Completely  Integrated 
Reference  Instrumentation  System  (CIRIS)  used  to  test  INS’s  at  Holloman  AFB,  NM 
[13:  1-1].  The  analysis  included  a  twenty-two  state  DGPS/INS  (GPS/INS  with 
differential  corrections)  Kalman  filter  which  was  practical  for  airborne 
implementations.  A  summary  of  the  twenty-two  state  filter's  achieved 
performance  is  presented  in  Table  1-1.  [  13:  Appendix  II. 

Note  that  Table  1-1  only  presents  results  for  INS  states  of  interest  to  the  kill 
zone  estimation  problem.  Also  note  that  while  the  reported  acceleration  errors  are 
not  directly  available  from  Negast's  thesis,  they  are  considered  to  be  obtainable 
from  typical  INS  accelerometers.  Because  the  kill  zone  estimation  problem  has 
been  defined  with  a  constant  altitude,  planar  geometry,  only  north  and  east 
velocity  and  acceleration  errors  are  required.  Because  the  volume  cross-section  is 
always  defined  in  9t(t)  and  0(t)  coordinates  relative  to  the  friendly  aircraft. 
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knowledge  of  precise  aircraft  position  is  irrelevant  to  the  kill  zone  estimation 
problem,  and  the  applicable  INS  position  errors  have  been  omitted  from  Table  1-1. 
Vertical  position  remains  only  as  an  aid  to  estimating  the  scenario  altitude  for 
missile  drag  considerations.  The  discussion  of  all  modelling  considerations  is 
deferred  to  Chapter  11. 


Table  l-l  DGPS-INS  Performance  Results 


INS  State 

RMS  Errors 

North  Velocity 

<  0.1  ft/sec 

Fast  Velocity 

<0.1  ft/sec 

North  Acceleration 

50  x  10"k  x  true  value 

East  Acceleration 

50  x  10"6  x  true  value 

Vertical  Position 

Oft 

Rockwell  International  is  currently  testing  the  concept  of  using  a  guided 
sphere  to  intercept  a  missile  and  achieve  kinetic  kill  [151.  The  sphere  uses  an 
onboard  sensor  and  small  lateral  thrusters  to  provide  proportional  navigation 
guidance  corrections  for  small  errors  in  the  estimated  intercept  point.  The  use  of  a 
"smart"  projectile  with  lateral  thrust  simplifies  the  required  missile  trajectory 
estimation  problem.  Furthermore,  Rockwell's  initial  testing  of  the  guided  sphere 
has  been  primarily  limited  to  space  scenarios,  where  the  missile  flight  path  is  more 
stable,  time-to-go  (T^)  is  not  a  critical  parameter,  and  the  missile  flies  a 

preprogrammed  ballistic  trajectory  (i.e.  the  missile  is  not  reacting  to  the  uncertain 
movements  of  an  evading  aircraft). 

This  thesis  differs  from  the  Rockwell  efforts,  in  that  this  thesis  assumes  the 
use  of  an  unguided  spread  of  projectiles  to  accomplish  the  missile  kill.  Successfully 
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using  a  "dumb"  projectile  necessarily  requires  accurate  estimation  of  the  future 
missile  trajectory.  It  has  been  suggested  that  a  priori  knowledge  of  future  aircraft 
maneuvers  is  key  to  investigating  the  kinetic  kill  zone  problem  (201.  The 
combination  of  special  aircraft  maneuvers,  along  with  accurate  estimation  of  the 
future  flight  trajectory  of  a  homing  missile,  may  enable  one  to  force  the  missile 
onto  a  path  which  optimizes  the  kinetic  kill  zone. 

1 .4  Assumptions 

This  thesis  assumes  that  the  encountered  homing  missile  uses  proportional 
navigation  or  a  derivative  thereof  (including  pure  pursuit,  lead-angle,  constant 
bearing  etc.)  as  the  primary  guidance  law,  regardless  of  the  type  of  seeker  head 
employed  (such  as  active/semi-active  radar,  infrared,  home-on-jam,  optical  etc.). 
This  thesis  exploits  the  known  attributes  of  proportional  navigation  guidance 
independent  of  the  region  of  the  electromagnetic  spectrum  in  which  the  missile 
seeker  head  operates. 

This  thesis  further  assumes  that  available  Missile  Warning  Receiver  (MWR) 
technologies  will  provide  accurate  initial  estimates  of  four  parameter  sets.  First, 
Tune-to-Go  (T-r^),  provides  the  time  remaining  before  missile  impact  with  the 
evading  aircraft.  The  TTC  provided  must  be  large  enough  to  allow  problem 
execution.  (What  defines  this  sufficient  time  is  an  output  parameter  of  this  thesis). 
Second,  the  MWR  verifies  existence  of  a  tracking  missile  which  starts  the  scenario. 
Third,  the  MWR  estimates  the  initial  angular  position  of  the  tracking  missile.  This 
is  particularly  important,  since  a  ladar  has  excellent  angular  resolution  --  less 
than  50  i*rad  has  been  demonstrated  —  but  is  very  poor  at  initial  acquisition  in  a 
large  field  of  view.  Fourth,  the  MWR  provides  initial  estimates  for  other  filter 
states  and  covariances  as  available.  Such  MWR  technology  is  currently  under 
development  by  Wright  Laboratory's  Avionics  Directorate's  Electronic  Warfare 
Division,  ECM  Advanced  Development  Branch,  Electro-optics  Group  (WL/AAWD-2). 
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Mr  Kevin  McKamey,  program  manager  for  WiyAAWD-2's  current  MWR 
development,  has  provided  the  initial  condition  values  summarized  in  Table  1-2  as 
being  valid  based  on  expected  MWR  technology  developments  [12]. 

Restrictions  are  enforced  on  the  allowed  maneuvers  of  the  evading  aircraft. 
Two  justifications  are  given  for  this  assumption.  The  missile  is  primarily  being 
controlled  by  its  proportional  navigation  guidance  law  which  commands  the 

Table  1-2  MWR  Provided  Initial  Conditions 


Initial  Parameter 

RMS  Error 

Range 

3.0  meters 

Range  Rate 

5.0  meter/ sec 

Angle 

2.0mrad 

Time-to-go 

10%  error 

missile  to  intercept  the  target  aircraft.  Therefore  the  aircraft  trajectory  becomes 
in  essence  a  series  of  guidance  and  control  commands  given  to  the  missile.  If  the 
aircraft  were  allowed  to  perform  a  random  maneuver  during  the  TTI,  the  aircraft 
maneuver  would  become  equivalent  to  a  random  guidance  command  given  to  the 
missile.  However  if  the  allowed  maneuvers  of  the  aircraft  are  benign  such  that 
they  can  be  accurately  estimated  by  an  extended  Kalman  filter  a  priori,  the 
maneuvers  instead  become  a  control  command  given  to  the  missile  to  force  a 
favorable  outcome. 

While  an  easily  predictable  maneuver  is  somewhat  restrictive,  the 
restrictions  are  considered  justifiable  given  that  this  thesis  represents  a  first 
approach  to  the  kill  zone  problem.  As  such  this  thesis  focuses  on  modeling  and 
analyzing  the  kill  zone  volume,  and  not  on  optimally  modeling  an  evading  aircraft. 
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Assuming  that  the  kill  /one  can  adequately  be  predicted  « lih  readily  predictable 
aircraft  maneuvers,  more  defensive  maneuvers  can  be  added  later  provided  they 
meet  the  following  condition:  The  filter  must  have  an  adequate  means  of 
consistently  and  accurately  estimating  apriori  the  aircrafts  maneuver.  Selecting  a 
benign  aircraft  trajectory  for  this  study  fulfills  that  condition.  As  a  final  note, 
restricting  the  allowed  aircraft  maneuvers  does  not  equate  to  exact  knowledge  of 
aircraft  maneuvers  actually  being  executed  (which  can  never  be  truly  known), 
but  rather  serves  to  enhance  the  accuracy  to  which  an  aircraft  maneuver  is 
known.  Sensitivity  of  the  final  predicted  kill  zone  volume  to  this  known  maneuver 
precision  is  analyzed. 

Since  this  thesis  depends  heavily  on  the  use  of  previous  thesis  studies,  most 
notably  Cusumano  and  Deponte,  and  Negast,  it  assumes  the  models  developed  in 
those  thesis  studies,  and  the  results  obtained  are  correct. 

There  is  currently  no  known  accepted  method  of  determining  a  homing 
missile's  T^  in  a  dynamic  environment  involving  evasive  aircraft  maneuvers. 

The  T-jc  problem  is  largely  driven  by  unknown  and  undeterminable  future 
maneuver  relations  between  the  missile  and  aircraft.  As  such,  Tj£  as  previously 
defined  is  always  changing.  In  order  for  a  predicted  volume  at  a  predicted  Tj-j  to 
make  physical  sense,  a  good  estimate  of  T^  needs  to  be  available.  For  purposes  of 
solving  the  thesis  problem,  an  overly  simplistic  TTC  model  is  developed  and  used 

with  the  understanding  that  such  technology  is  under  development  [12].  This 
thesis  treats  TjG  as  an  analysis  parameter  so  that  the  results  presented  are 
independent  of  current  T-jg  estimation  problems. 

1.5  Scope 

This  thesis  begins  with  the  detection  of  a  homing  missile  already  in  flight 
that  is  tracking  the  friendly  aircraft.  Initial  conditions  in  Table  1-2  are  provided 
by  the  MWR  and  include  a  valid  missile  detection,  angular  direction  (in  milli- 
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radians),  Tr-,  and  initial  filter  conditions  for  all  state  values  an  MWR  would  require 
to  process  detection,  position,  and  T^;.  The  thesis  model  predicts  an  intercept  point 
and  a  bounded  volume  at  a  specified  remaining  and  T  p|.  The  central  intercept 

location  of  the  kill  zone  v  olume  is  in  general  not  a  single  point  but  rather  a 
solution  which  generates  the  point  based  on  key  input  parameters  such  as  T-ip,, 
initial  position,  final  Tpi  ,  allowed  aircraft  maneuvers,  and  measurement  rates,  as 
well  as  all  the  associated  error  tolerances  for  these  parameters.  This  thesis  does  not 
address  any  of  the  kinematic  problems  associated  with  the  actual  deployment  of  the 
kinetic  kill  cannister.  For  purposes  of  this  thesis,  a  prediction  of  a  "good  volume 
cross-section"  is  assumed  as  lethal  to  a  homing  missile.  This  thesis  makes  maximum 
use  of  existing  missile  trajectory,  internal  missile  dynamics,  and  INS  error  models, 
and  Kalrr  an  filter  models  developed  in  previous  thesis  studies.  This  thesis  focuses 
on  analyzing  the  kinetic  kill  zone  volume  and  determines  how  changing  each  of 
the  assumed  input  parameters  changes  the  output  volume  parameters. 

Figure  1-4  shows  a  block  diagram  representation  of  the  thesis  model.  The 
dynamics  equations  describing  the  desired  aircraft  maneuvers  (a)  are  programmed 
into  the  Kalman  filter's  state  equations.  Pilot  deviation  from  the  filter- 
programmed  trajectory,  is  represented  by  an  appropriate  noise  strength  driving 
the  extended  Kalman  filter's  aircraft  state  covariances.  The  true  aircraft  trajectory 

(A)  is  generated  assuming  a  constant  speed/constant  angle-rate  turn  driven  by 
pilot-induced  noise.  Jinking  maneuvers  are  not  used  since  they  would  tend  to  be 
unpredictable  by  the  filter  and  would  therefore  degrade  the  final  solution.  Noise 
corrupted  measurements  of  the  true  aircraft  trajectory  (b)  are  provided  to  the 
filter  from  the  INS.  The  true  missile  dynamics  (C,D)  are  simulated  using  a  modified 
version  of  Cusumano  and  Deponte's  5th  order  missile  model  [5].  The  ladar  takes 
measurements  between  the  true  missile  position  (E)  and  the  true  aircraft  position 

(B)  to  provide  the  extended  Kalman  filter  with  noise  corrupted  measurements  of 


1-13 


. .  . .  —  - 

i 

range,  range  rate,  angle,  and  angle  rate  (f).  The  filter  then  updates  current  filter 
state  estimates  with  the  ladar  and  INS  measurements  (g)  and  propagates  the  states 
forward  one  measurement  period  to  generate  new  state  estimates  for  the  next 
measurement  cycle.  An  indcpth  analysis  of  these  system  functions  is  provided  in 
Chapter  11. 

1.6  Approach  and  Methodology 

The  Kalman  filter  is  simulated  using  a  Multimode  Simulation  for  Optimal 
Filter  Evaluation  (MSOFF)  software  [31.  An  extended  Kalman  filter  model  is  used  to 
estimate  the  nonlinear  truth  model  states.  The  system  truth  model  generates 
aircraft  flight  dynamics  internally.  The  simulation  assumes  filter  access  (during 
measurement  cycle)  to  a  GPS-aided  INS  based  on  a  twenty-two  state  Kalman  filter 
comparable  to  that  presented  in  the  Negast  thesis  [131.  All  computer  simulations 
are  run  using  Standard  FORTRAN  77  code. 

There  are  four  major  parts  of  the  extended  Kalman  Filter  model:  ( 1 )  the 
friendly  aircraft  INS  states,  (2)  the  missile  flight  dynamics  model,  (3)  the  ladar 
measurement  model,  and  (4)  the  missile  body  dynamics  model.  The  final  volume 
cross-section  is  generated  by  propagating  the  filter's  steady  state  kinematics 
estimation  results  forward  by  the  T-n  parameter. 

The  thesis  uses  a  parametric  approach  to  analyzing  the  various  elements. 
First  a  general  baseline  performance  model  is  established  using  low  (or  zero  where 
appropriate)  noise  strengths  and  simple  aircraft  and  missile  trajectories.  Specific 
problem  parameters  of  interest  are  varied  from  the  baseline  model  to  determine 
kinetic  kill  zone  sensitivity  to  that  parameter.  Then  multiple  parametric  studies 
are  accomplished  using  more  complex  scenarios,  to  determine  how  the  parameters 
of  interest  inter-relate.  Parameters  to  be  analyzed  include  noise  strengths  used, 
aircraft  trajectories,  missile  approach  angles,  required  Tjj,  and  missile  dynamics. 
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Fig  1  -4  Overview  of  Simulation  model 


1.7  Summary 

Chapter  l  introduced  and  defined  the  kinetic  kill  zone  problem  under  study 
in  this  thesis.  Section  1 . 1  presented  the  background  rationale  lor  performing  this 
study,  with  a  specific  problem  statement  and  established  research  goals  following 
in  Section  1.2.  Previous  research  and  other  related  work  which  made  significant 
contributions  to  this  study  was  summarized  in  Section  1.3.  Section  1.4  presented  the 
basic  assumptions  underlying  the  problem  definition  as  presented  in  Section  1.2, 
and  gave  rationale  supporting  those  assumptions  where  applicable.  Section  1.5 
scoped  the  work  to  be  accomplished  in  this  thesis  study  and  overviewed  the 
simulation  model  developed  in  detail  in  Chapter  II.  Section  1.6  concluded  the 
chapter  with  a  discussion  of  simulation  methods,  and  parametric  analysis 
techniques  used  to  accomplish  the  established  research  goals. 

Chapter  II  develops  mathematical  equations  for  the  filter  and  truth  models 
used  in  this  study.  System  and  filter  states  are  established.  System  dynamic  and 
measurement  noise  models  are  developed.  Finally  any  special  initial  conditions 
(mean  or  covariance)  used  by  the  models  are  presented  with  supporting 
justification. 

Chapter  III  overviews  the  simulation  runs  conducted  to  parametrically 
analyze  the  kinetic  kill  zone  estimation  problem.  The  simulations  focus  on  three 
primary  areas;  the  missile-aircraft  scenario  geometry,  variation  of  the  assumed 
ladar  parameters,  and  variation  of  the  assumed  aircraft  parameters.  High  and  low 
performance  baselines  are  established  for  the  aircraft  and  ladar  parametric 
studies. 

Chapter  IV  presents  the  simulation  results.  Particular  attention  is  focused 
on  how  the  predicted  volume  varies  as  a  function  of  range,  time,  and  the 
parameters  under  study.  Cross-section  figures  of  merit  are  established  as  a  means 
of  standardizing  the  analysis.  Best  case  and  worst  case  baselines  are  established  for 
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performance  comparisons.  Parameters  which  had  either  very  large  or  negligible 
influence  on  the  predicted  volume  are  presented. 

Chapter  V  presents  the  thesis  study's  primary  conclusions.  The  kinetic  kill 
zone  prediction  approach  presented  in  this  thesis  (i.e.  extended  Kalman  filtering 
using  propagated  covariance  estimates  to  predict  the  cross-section  statistics)  is 
analyzed  in  terms  of  its  feasibilty  based  on  the  Chapter  IV  results. 
Recommendations  for  continued  research  are  also  presented  here. 
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II  The  Simulation  Models 


2.0  Introduction 

Chapter  II  develops  the  models  presented  in  Figure  1-4  in  detail.  Discussion 
starts  with  the  missile  and  aircraft  dynamics.  Appropriate  coordinate  frames  are 
established  for  the  missile  and  aircraft  dynamics  in  Section  2. 1 .  Kinematic 
equations  of  motion  are  then  developed  using  the  established  coordinate  frames  in 
Section  2.2.  Special  considerations  such  as  missile  drag,  internal  missile  dynamics, 
and  aircraft  maneuver  dynamics  are  also  covered  in  Section  2.2.  Section  2.3  follows 
with  a  discussion  of  truth  states  and  filter  states  selected  from  the  dynamics 
equations  developed  in  Section  2.2.  Section  2.4  discusses  the  noise  terms  used  in  the 
truth  and  filter  models,  including  system  process  noise,  filter  Q.  (dynamic  noise 
strength)  values,  sources  of  ladar  measurement  noise,  and  INS  measurement  noise. 
(Note  that,  in  this  context,  INS  measurement  noise  is  taken  to  mean  the  DGPS/INS 
filter's  uncertainty  about  state  values  when  the  INS  filter  is  requested  to  supply  an 
INS  update  to  the  tracking  filter.  Since  an  INS  filter  model  was  not  developed  for 
this  thesis,  the  tracker-requested  INS  updates  are  processed  as  measurements  from 
an  external  source.)  The  chapter  concludes  with  a  discussion  of  initial  conditions 
presented  in  Sections  2-4  and  2-5. 

2.1  Reference  Frames 

There  are  four  reference  frames  used  in  developing  the  system  dynamics 
model.  The  four  reference  frames  illustrated  in  Figure  2-1  (a-b)  include  the 
aircraft's  and  missile's  reference  frame  (a),  the  aircraft- to-missile  relative 
coordinate  frame  (b),  and  the  missile-to-aircraft  relative  coordinate  frame  (b).  The 
aircraft  receives  friendly  aircraft  kinematic  measurement  updates,  defined  in 
north-east-down  (NED)  coordinates,  (the  aircraft  reference  frame)  which  are 
supplied  by  an  onboard  Inertial  Navigation  Subsystem  aided  with  Global 
Positioning  Satellite  updates  and  Differential  Corrections  (INS/DGPS).  The 
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aircraft-to-missile  reference  frame  is  defined  in  (Rt,0t)  coordinates  originating  at 
the  aircraft.  Relative  missile-io-aircrali  kinematic  measurements  are  supplied  by 
the  ladar  to  the  aircraft's  Kalman  filter  in  this  reference  frame.  The  missile-to- 
aircraft  reference  frame  is  defined  in  (Rm  =  Rt  ,0m=0t  +  *)  coordinates.  The 
missile  uses  relative  kinematic  measurements  obtained  in  the  missile-to-aircraft 
coordinate  system  to  generate  the  appropriate  lateral  accelerations,  Aj.,  which 
guide  the  missile  to  the  aircraft.  Since  no  direct  access  is  available  to  the  missile's 
kinematic  states,  the  missile's  reference  frame  is  free  to  be  chosen  in  any 
appropriate  coordinates  and  the  NF.D  coordinates  used  by  the  aircraft  are  chosen 
for  convenience. 

Since  this  thesis  study  is  limited  to  an  in-the-plane  analysis,  all  references 
to  NED  coordinates  are  assumed  to  be  at  constant  altitude  (down  reference  is  fixed 
constant)  and  all  references  to  (R,0)  coordinates  are  assumed  to  lie  in  a  constant 
altitude  plane. 

2.2  Developing  The  Simulation  Model 

Development  starts  with  an  aircraft  flying  a  specified  trajectory  in  the 
aircraft  coordinate  frame.  A  missile  is  assumed  to  be  tracking  the  aircraft  with 
sufficient  closing  velocity  to  achieve  intercept.  The  missile  is  positioned  in  (Rm 
0)n)  coordinates,  using  the  missile-to-aircraft  reference  frame.  0mis  defined  as 
being  counterclockwise  from  east  in  the  aircraft  reference  system.  The  respective 
aircraft  and  missile  headings  are  defined  as 

Rm=ARCTAN(Vmn/VItE)  (2-1) 

6t  =  ARCTAN  ( Vtn  /  Vte)  ( 2-2 ) 


where 

ftm  is  the  missile  heading, 

&t  is  the  aircraft  heading, 

Vmn  is  the  missile's  north  velocity  component, 
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Vfljg  is  the  missile's  east  velocity  component, 

Vtn  is  the  aircraft’s  north  velocity  component, 

Vte  is  the  aircraft's  east  velocity  component. 

Since  the  missile  reference  frame  is  coupled  to  the  aircraft  reference  frame 
by  the  commanded  lateral  accelerations,  acceleration  must  necessarily  be 
propagated  as  a  state.  To  this  effect,  equations  of  acceleration  are  taken  from 
Symon  which  describe  acceleration  in  (91,0)  coordinates  [17]. 

ar  =  d29t/dt2  -  91(d0/3t)2  (2-3) 

a0  =  9l(d20/3t2)  +  2  (391/at)  (30/at)  (2-4) 

where 

ar  is  the  radial  acceleration, 
a0  is  the  tangential  acceleration. 

Using  the  state  definitions  provided  later  in  Table  2-2,  (for  a  discussion  on  the 
selection  of  filter  and  system  states  see  Section  2.3)  a  coordinate  transformation 
between  (X,Y)  -or-  NED  and  (91,0)  coordinates,  and  properties  of  vectors,  the  missile 
and  aircraft  NED  referenced  accelerations  are  rewritten  in  the  missile  to  aircraft 
reference  frame,  using  the  coordinate  transformation  matrix. 

rx'j  =  r  cosp  sin  p  i  r  xi  (2-5) 

lY'J  1  -  SIN  p  COSp  J  lY  J 

where  p  =  0m  -  (ftm  or  ftt  respectively), 
and 

fcm  is  the  missile  heading, 
ftt  is  the  aircraft  heading, 

0m  is  the  missile  to  aircraft  line-of-sight. 

The  transformed  missile  acceleration  (from  NED  to  (91,0)  coordinates)  becomes 

laj  =  ( (aV^/at)2  +  (aVjng/at)2  )1/2  (2-6a) 
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arm=  lamICOS(  0m  -  ARCTAN(  (aVmn/at)/(  V^/at)))  (2- 6b) 

a0m=  -  laJSIN(  0m-  ARCTAN(  (.W^/rHl/UiV^/at)))  (2-6c) 

where 

la^  is  the  magnitude  of  the  missile  acceleration  vector, 

a,.  is  the  radial  missile  acceleration, 
rm 

a^is  the  tangential  missile  acceleration. 

Using  the  same  transformation  on  the  aircraft  acceleration  vector  yields 

latl  =  ( (aVtn/at)2  +  (dVte/at)2  )l/2  (2-7a) 

art  =  latICOS(  0m  -  ARCTAN1  OVtn/at)/(dVte/at))  )  (2-7b) 

a0{  =  -  latISIN(  0m  -  ARCTAN(  (3Vtn/at)/(aVte/dt))  )  (2- 7c) 


where 

latl  is  the  magnitude  of  the  aircraft  acceleration  vector, 
art  is  the  radial  aircraft  acceleration, 
a^  is  the  tangential  aircraft  acceleration. 

The  relative  missile-to-aircraft  acceleration  vector  in  the  missile-to-aircraft 
reference  frame  is  derived  by  subtracting  Equation  set  (2-6)  from  Equation  set  (2- 
7)  to  yield: 


«ir  *.  — 

rmt  rm  U 
a0mt =  a0m '  a0t 


(2-8) 

(2-9) 


Combining  Equations  (2-8)  and  (2-9)  with  Equations  (2-1)  and  (2-2)  yields  the 
following  expressions  for  radial  acceleration  and  angular  acceleration: 


d2m/dt2  =  ar  -  a_  +m(d0/d t)2 
rt  rm 

(2-10) 

a20/at2  =  (a^  -a0m  -2(a3ft/at)(d0/at))/JR 

(2-11) 

This  concludes  the  development  of  the  kinematic  model.  The  section 
continues  with  a  development  of  missile  and  aircraft  dynamics.  State  selection  is 
deferred  to  Section  2.3. 
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2.2.1  The  5th  Order  Missile  Model 

Cusumano  and  De  Ponte  developed  a  5th  order  missile  truth  model  which  was 
used  as  a  performance  benchmark  against  which  to  evaluate  lower  order  filter 
models  shown  in  Figure  2-2  [5|.  The  truth  model  autopilot  had  time-varying  poles 
and  zeros  which  altered  the  overall  autopilot  transfer  function  slightly  with  time. 
The  5*  order  missile  model  was  selected  to  represent  the  truth  model's  internal 
missile  dynamics  for  this  thesis.  However,  the  time-varying  poles  and  zeros  were 
replaced  by  constant  values  equal  to  the  time- zero  values.  This  was  done  based  on 
the  recommendation  of  Cusumano  and  De  Ponte,  since  the  time-varying  transfer 
function's  dominant  poles  and  zeros  did  not  change  significantly  from  their  time- 
zero  values,  and  the  time-zero  autopilot  transfer  function  was  always  a  good 
approximation  to  the  time  varying  autopilot  transfer  function  [5]. 

2.2.2  Series  To  Parallel  Conversion 

The  autopilot  portion  of  the  5th  order  missile  model  is  implemented  as  a 
decoupled  transfer  function,  shown  in  Figure  2-3,  by  doing  a  partial  fraction 
expansion  on  the  autopilot  transfer  function  [7].  The  autopilot  was  converted  in 
order  to  decouple  the  autopilot  poles  and  zeros  in  the  MSOFE  computer  code. 
Decoupling  the  states  allowed  pole  and  zero  values  to  be  altered  more  easily.  The 
residues  of  the  decoupled  states,  listed  in  Table  2-1,  are  evaluated  as 

r,  =(s  +  pj)  H(s)  |s  =  -Pi  (2-12) 

where 

rt  is  the  residue  of  pole  pp 

Pj  is  the  current  pole  under  evaluation, 

H(s)  is  the  5th  order  missile’s  transfer  function, 
s  is  a  Laplace  transform  variable. 


2-6 


3 


' _ 

cn 

co 

prj 

vO 

00 

oo 

+ 

ro 

f— H 

cn 

O 

r—4 

H 

V 

on 

rsi 

H 

P 

P 

ll 

ll 

if) 

\ 

\ 

ro 

i—4 

i-h 

a 

a 

i 


2-; 


ps  =  63.34 
zi  =  33.34 

Z2  =31.10  Figure  2-2  5th  Order  Cusumano  and  Deponte  Missile  Guidance  Strip 
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Figure  2-3  Paraliel  Implementation  of  5th 
Order  Missile  Transfer  Function 


Table  2- 1  Converted  Autopilot  Poles 


5th  Order  Missile 

Pole 

Converted  Residue 

Values 

s  +  3 .33 

3.33  (no  change) 

s  +  13.33 

13.33  (no  change) 

s  +  3.06 

-.098751  (r3) 

s  + 21.88 

+  .783663  (r4) 

s  +  b3.34 

+  1.12986  (r5) 

DC  Gain 

-3.964 

2.2.3  The  Missile  Guidance  Law 

The  missile  generates  lateral  acceleration  (AjJ  in  response  to  aircraft 
maneuvers  according  to  the  proportional  navigation  guidance  law 

ALcomd(s>  =  n  *»(s)  s0(s)  (2-13a) 

AL^®)  =  ALcomd(s)  H(s)  (2-13b) 

where 

Al(s)  is  the  missile's  generated  lateral  acceleration, 

ALcomd(s)  is  the  guidance  law's  commanded  lateral  acceleration, 
n  is  a  constant  of  proportional  navigation  between  3  and  6  [22], 
s9t(s)  is  the  range-rate  between  the  missile  and  aircraft, 
s0(s)  is  the  line-of-sight  rate  between  the  missile  and  aircraft, 

H(s)  is  the  5lh  order  missile  transfer  function. 

2.2.4  Atmospheric  Effects 

The  missile  used  in  this  simulation  is  assumed  to  be  in  a  post-burnout,  coast- 
to-intercept  phase  of  flight  The  missile  is  therefore  assumed  to  be  constant  mass 
and  kinetically  limited.  The  missile  has  two  drag  terms  which  significantly  affect 
the  missile's  future  flight  trajectory.  Both  the  coefficient  of  drag,  CD  and  the 
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coefficient  of  lift,  C^,  generate  drag  force  acting  in  opposition  to  the  missile's 
current  velocity  vector.  However,  the  coefficient  of  lift  also  produces  a  lift  force 
perpendicular  to  the  missile  velocity  vector.  The  amount  of  lift  force  generated  is 
controlled  by  the  missile's  guidance  law  and  transfer  function  from  Equation  (2- 
13).  The  amount  of  additional  forward  drag  induced  by  the  coefficient  of  lift  is  a 
function  of  the  lift  force  requested. 

The  detailed  development  of  the  two  drag  coefficients  is  covered  thoroughly 
in  Zarchan  [22],  but  is  summarized  here  for  completeness.  The  coefficient  of  drag 
is  a  function  of  four  deterministic  parameters  -  the  air  density  (p),  the  missile 

surface  area  (cross-section)  to  mass  ratio  (S/M),  the  speed  of  sound  (y),  the 
scenario  altitude  (h)  -  and  the  coefficient  of  zero  lift  (Cqq).  The  coefficient  of  lift 

is  a  function  of  the  air  density  (p),  the  missile  surface-to-mass  ratio  (S/M),  the 
missile  velocity  (Vm)  and  the  commanded  lateral  acceleration  (AIromci)  [22].  The 

drag  equations  are 

air  density 
coefficient  of  drag 
coefficient  of  lift 
drag  acceleration 

where 

ALcomd  is  commanded  lateral  acceleration, 

A[)(t)  is  the  total  missile  drag  due  to  lift  induced  and  normal  drag, 

Qx)  is  the  no-lift  induced  drag, 

Cql  *s  the  lift  induced  drag. 

Vm  is  the  missile  velocity  (magnitude). 


P  =  0.002378  e-(h/30’000>  (slg/ft3) 

(2-14) 

Cq=  03  p  S/M  Cqo 

(2-15) 

Cl  =  2  Alcoi,1ci/(S/M  p  Vm^) 

(2-16) 

AD(t)  =  -(  CD+CDL)Vm2 

(2-17) 
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Cusumano  and  De  Ponte  further  develop  the  and  terms  in  their 
thesis  work  [51: 

zero  lift  drag  (’lX)  =  2  y  1/2  /  Vm1/2  (2-18) 

lift  induced  drag  (^L  =  4ALcomd2/(p2  (S/M)2  Vm^  Y)  (2-19) 

where  y  is  the  speed  of  sound  in  air. 

2.2.5  Aircraft  Maneuver  Generator 

All  aircraft  maneuvers  are  assumed  to  be  benign  for  this  simulation  for  two 
reasons.  When  the  filter  switches  from  the  track  mode,  in  which  the  filter 
continues  to  converge  the  state  covariances,  to  the  predict  mode,  the  filter 
propagates  the  state  covariances  forward  by  the  desired  intercept  time  (Tjj). 

Because  the  filter  accesses  a  priori  information,  no  measurements  from  either  the 
ladar  or  the  INS  are  accessible  during  this  extended  cycle.  All  anticipated  aircraft 
maneuvers  act  as  control  inputs  to  the  missile's  flight  trajectory,  as  coupled 
through  the  proportional  navigation  guidance  laws.  The  filter  estimate  of  the 
aircraft  trajectory  therefore  accounts  for  two  processes  through  white  noise  of 
strength  Oft).  The  noise  strength  reflects  limitations  in  the  filter's  aircraft 
trajectory  model,  and  the  pilot's  limited  ability  to  execute  the  desired  maneuver 
without  error.  A  benign  trajectory  simplifies  the  analysis  and  allows  the  Inter  to 
model  the  nominal  aircraft  trajectory  accurately.  The  benign  trajectory  is  also 
easier  for  the  pilot  to  track,  resulting  in  an  appropriate  decrease  in  Oft). 

The  filter  equations  are  developed  as  follows  assuming  a  constant  speed 
constant  angle-rate  turn  [17]: 


1  acceleration! 

Vf2/Rtum  =  9’8(3ievei  -  u>2  Rtum 

(2-20) 

position 

Xte(t)  =  Rtum  GOS(ot+  8) 

(2-2  la) 

Ytn(t)  =  Rtum  SIN(o)t+  8) 

(2-2  lb) 

velocity 

Vte(t)  =  -o)RturnSIN(cot  +  8) 

(2-22a) 

Vtn(t)  =  aiRturnCOS(cBt  +  8) 

(2-22b) 
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acceleration 


acceleration 

Atf,(t)  =  -oi^RmrnGOSOot  +  0) 

(2-2  3a) 

A,n(l)  =  -u)/RlurnSIN(o>t  +  0) 

(2-22b) 

jerk 

Iteft)  =  <'»-;*RturnSIN(",l+  0) 

(2-24a) 

Jtn(l)  =  -«.,3RturnCOS(<»l  +  0) 

(2-24b) 

where 

Vt  is  the  aircraft's  velocity  (magnitude), 

Gievel  *s  ttiH1  acceleration  in  g's, 

RtUm  ts  the  turn  radius, 

oi  is  the  angle-rate  of  the  turn  in  rad/ sec, 

8  is  a  phase  angle  determined  by  the  aircraft's  initial  velocity  heading, 
and  all  kinematic  terms  are  derived  in  the  NED  reference  frame. 

The  equations  of  circular  motion  (Equations  (2-20)  through  (2-24))  are  selected  to 

model  the  aircraft  motion  for  the  following  reasons:  ( 1 )  The  trajectory  is  easily 

modeled  mathematically  by  both  the  filter  and  the  truth  system;  (2)  The  trajectory 

is  easily  executed  by  a  pilot;  (3)  The  trajectory  meets  qualifications  of  being  a 

benign  "easily  predictable"  trajectory;  (4)  The  trajectory  is  a  subset  of  potential 

evasive  maneuvers  which  may  be  included  as  part  of  a  follow-on  analysis. 

2.3  State  Selection 

The  state  selection  process  ( for  both  the  truth  model  and  the  filter  model) 
begins  with  a  summary  of  the  three  key  facts  which  dictated  the  choice  of 
reference  frames.  1 )  All  aircraft  INS  measurements  occur  in  the  NED  reference 
frame.  2)  All  aircraft  ladar  measurements  occur  in  the  aircraft-to-missile 
reference  frame  (i.e.  polar  coordinates  originating  from  the  aircraft).  3)  The 
missile  navigates  using  target  information  obtained  in  the  missile-to-aircraft 
reference  frame  (i.e.  polar  coordinates  originating  at  the  missile).  Because  the 
most  complicated  dynamics  occur  in  the  missile-to-aircraft  reference  frame  (i.e. 
the  missile  guidance  commands),  it  is  chosen  as  the  primary  reference  frame  for 
the  truth  and  filter  models. 
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2.3.1  The  Truth  States 

The  state  selection  process  begins  with  Equations  (2-10)  and  (2-11),  which 
describe  radial  and  angular  acceleration,  respectively,  in  the  missile-to-aircraf't 
reference  frame.  From  the  two  equations  ten  states  are  selected.  They  include 
range,  range-rate,  angle,  angle-rate,  north  missile  velocity,  east  missile  velocity, 
north  aircraft  velocity,  east  aircraft  velocity,  north  aircraft  acceleration  and  east 
aircraft  acceleration.  The  four  aircraft  states  are  measurable  directly  from  the 
aircraft  INS.  The  ladar  provides  relative  aircraft-to-missile  measurements  of 
range,  range-rate,  angle  and  angle-rate.  The  missile  velocity  states  are  measured 
by  combining  the  absolute  aircraft  measurements,  available  from  the  INS,  with  the 
relative  velocity  measurements  available  from  the  ladar.  All  ten  kinematic  states 
are  therefore  observable  with  respect  to  the  assumed  INS  and  ladar  measurements. 
Missile  velocity  is  propagated  as  acceleration  driven  by  drag,  Equations  (2-18)  and 
(2-19),  and  (guidance  law  commanded)  lateral  acceleration.  Equation  (2-13).  As 
seen  from  Equations  (2-18)  and  (2-19),  drag  is  a  function  of  absolute  missile 
velocity,  hence  the  reason  for  propagating  both  absolute  and  relative  missile 
kinematics. 

The  missile  transfer  function,  H(s),  contains  five  internal  states  (ref.  Figure 
2-3)  which  are  selected.  The  proportional  navigation  constant,  n,  is  a  parameter. 
From  Equations  (2-18)  and  (2-19),  two  additional  parameters  are  chosen  -  the 
scenario  altitude,  h,  and  the  missile's  surface  to  mass  ratio  (S/M).  No  additional 
states  are  required  to  model  missile  drag. 

Aircraft  kinematics  are  taken  from  Equations  (2-20)  thru  H  4b).  Two 
new  states  representing  aircraft  velocity  are  added  along  with  g-level  which  is  a 
parameter.  RtUm>  (0>  and  8  can  be  calculated  internally  from  the  states  and 

parameters  already  available  and  are  not  needed  as  separate  states.  A  summary  of 
selected  truth  model  states  is  provided  in  Table  2-2A. 
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Table  2-2A  Truth  Model  States 


TRUTH 

STATE 

DEFINITION  OF  STATE  VARIABLE: 

UNITS  OF  STATE 
VARIABLE 

True  Missile  Velocity  X  Direction 

m  /  sec 

K6H 

True  Missile  Velocity  Y  Direction 

m  /  sec 

Xs(3) 

True  Line-of-Sight  Angle 

radians  (rad) 

Xs(4) 

True  Line-of-Sight  Rate 

rad  /  sec 

Xs(5) 

True  Missile-to-Aircraft  Range 

meters  (m) 

Xs(G) 

True  Range-Rate 

m  /  sec 

Xs(7) 

Aircraft  East  Velocity 

m  /  sec 

Xjj(8) 

Aircraft  North  Velocity 

m  /  sec 

X*(9) 

Aircraft  East  Acceleration 

m  /  sec2 

Xs(10) 

Aircraft  North  Acceleration 

m  /  sec2 

mm 

Seeker  Bandwidth 

rad  /  sec 

Xs(12) 

Seeker-Commanded  Lateral  Acceleration 

m  /  sec2 

Xs(13)-(15) 

Missile  Autopilot  Transfer  Function 

m  /  sec2 

Al 

Generated  Lateral  Acceleration 

2Xs(13)  ~Xs(15) 

n 

Pro-Nav  Parameter  of  Missile 

const  (Val  =  3.0  -  6.0) 

h 

Scenario  Altitude 

feet 

S/M 

Missile's  Surface-to-Mass  Ratio 

m2  /  kg 

2.3.2  The  Filter  States 

Filter  states  are  chosen  identical  to  the  truth  states  except  as  noted  here.  All 
variables  listed  as  parameters  (except  g-level)  in  the  truth  model  are  considered  as 
unknown  constants  by  the  filter  and  are  therefore  modeled  as  the  outputs  of 
integrators  driven  by  white  pseudo-noise.  The  unknown  constants  include  n,  S/M, 
and  h.  In  addition,  the  unknown  missile  transfer  function  is  modelled  as  a  perfect 
direct  feedthrough  with  an  appropriate  increase  in  the  filter's  missile  dynamics 
noise  stmegth  and  in  the  pseudo-noise  strength  driving  n.  Although  the 
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simulations  are  set  up  to  be  strictly  in-the-piane  (contant  altitude)  scenarios, 
missile  drag  is  a  function  of  altitude  (albeit  a  weak  one).  It  is  not  expected  that 
altitude  estimations  will  have  a  significant  impact  on  filter  performance.  Altitude 
has  been  included  anyway  since  it  will  not  degrade  filter  performance.  If  one 
seeks  to  optimize  the  filter  design,  altitiude  estimation  can  be  removed  without 
significantly  changing  filter  performance.  The  selection  of  filter  states  is 
summarized  in  Table  2-2B. 

Table  2-2B  Filter  Model  States 


FILTER 

STATE 

DEFINITION  OF  STATE  VARIABLE 

UNITS  OF  STATE 
VARIABLE 

Xf(l) 

Est.  Missile  Velocity  X  Direction 

m  /  sec 

KH 

Est.  Missile  Velocity  Y  Direction 

m  /  sec 

Xf(3) 

Est.  Une-of-Sight  Angle 

radians  (rad) 

UEsm 

Est.  Line-of-Sight  Rate 

rad  /  sec 

m m 

Est.  Missile-to-Aircraft  Range 

meters  (m) 

Est.  Range-Rate 

m  /  sec 

Xf(7) 

Time-to-Go  Estimate 

seconds  (sec) 

mam 

Est.  of  ProNav  Parameter 

constant  4.5  nominal 

Xf(9) 

Est.  Surface-to-  Mass  Ratio 

m2  /  kg 

XflO) 

Est.  Missile  Lateral  Acceleration 

m  /  sec2 

warn 

Est.  Scenario  Altitude 

feet 

Xf(12) 

Aircraft  East  Velocity 

m  /  sec 

Xf(13) 

Aircraft  East  Acceleration 

m  /  sec2 

RSn 

Aircraft  North  Velocity 

m  /  sec 

Xf(15) 

Aircraft  East  Acceleration 

m  /  sec2 
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2.4  The  Noise  Terms 

All  noise  used  to  drive  the  model  ran  be  categorized  in  one  of  four  t>pes: 
system  process  noise,  INS  related  measurement  noise,  ladar  measurement  noise,  and 
filter  parameter  estimation  (or  tuning)  noise.  Process  noise  is  a  representation  of 
unmodelled,  or  partially  modelled  real  dynamic  effects.  Process  noise  always  refers 
to  noise  terms  driving  the  truth  model  states  Xs(t).  Tuning  noise  represents 

uncertainty  in  the  filter's  dynamic  model  and  can,  but  does  not  necessarily, 
represent  real  dynamic  effects.  ladar  measurement  noise  represents  limitations  in 
the  measurement  accuracy.  INS  measurement  noise  represents  the  INS  filter’s 
estimated  uncertainty,  which  is  propagated  in  the  INS  filter  as  a  state  covariance. 
The  kill  zone  estimation  filter  periodically  accesses  the  INS  filter,  to  update  the  kill 
zone  estimator's  propagated  aircraft  states.  The  INS  filter's  estimated  covariance 
becomes  the  kill  zone  estimation  filter's  measurement  noise  at  the  prescribed 
update  time. 

2.4.1  Process  Noise 

Process  noise  enters  the  model  through  the  truth  model  dynamic  equations 

as  [9] 

ax(t)/dt=JIt,s(tH  +  CHt,s(tHw(t)  (2-25) 

where 

X(t)  is  the  vector  of  states  previously  defined  in  Table  2-1, 

f{t,x(t){  is  a  vector  of  nonlinear  dynamic  equations  developed  in  Section  2-1, 

£Ht,x(t)}  is  a  vector  of  noise  shaping  terms, 

w(t)  is  a  vector  of  "white  noises"  which  partially  drive  the  system  dynamics. 
Aircraft  Process  Noise  The  aircraft's  acceleration  states  (  Xsgft),  Xsioft)  )  are 
given  a  process  noise  which  models  the  pilot's  inability  to  execute  a  perfect 
maneuver.  When  the  filter  switches  from  the  track  mode  (in  which  the  filter 
converges  the  covariance  estimates  for  each  state  to  its  steady  state  value)  to  the 
predictor  mode  (in  which  the  filter  tries  to  predict  the  missile  trajectory),  the  filter 
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loses  all  access  to  measurement  updates.  The  filter  assumes  a  nominal  aircraft 
trajectory  which  is  to  be  perfectly  executed.  Since  the  aircraft  trajectory 
represents  a  control  input  to  the  tracking  missile  (as  coupled  through  the  lateral 
acceleration  guidance  commands),  any  deviation  from  the  filter's  assumed 
trajectory  will  cause  a  corresponding  error  in  the  filter's  estimate  of  the  missile 
trajectory.  Most  pilots  are  capable  of  sustaining  constant  velocity,  constant  gee- 
level  maneuvers  within  ±0.1  g  for  at  least  30  second  intervals  [4],  a  much  larger 
time  interval  than  a  typical  kinetic  kill  zone  scenario  would  last.  Since  the  0. 1  g 
acceleration  error  is  controlled  and  bounded  by  the  pilot,  the  acceleration  error  is 
more  typical  of  a  measurement  noise  than  a  dynamic  process  noise.  A  noise 
strength,  Qft),  of  (0.1  gj^/sec  is  considered  an  adequate  approximation  to  represent 
the  pilot  induced  trajectory  error. 

Missile  Process  Noise  Figure  2-4a  shows  a  missile  seekerhead.  The  missile 
body  axis  is  shown  with  vector  (a),  the  missile  seeker  axis  with  vector  (b)  and  the 
target  LOS  with  vector  (c).  The  angle  between  vectors  (b)  and  (c)  is  the  missile’s 
tracking  error,  and  is  measured  by  the  missile  seeker.  The  tracking  error  is 
represented  on  the  horizontal  axis  of  Figure  2.5(a  &  b).  The  moving  target's  angle- 
rate  in  the  missile-to-aircraft  reference  frame  causes  a  change  in  vector  (b) 
relative  to  (c)  and  hence  a  change  in  the  tracking  error.  The  tracking  error  is 
passed  to  an  amplifier  (Figure  2.5)  which  generates  an  appropriate  angle-rate 
signal  to  null  the  tracking  error,  thus  allowing  the  missile  to  track  the  target.  The 
generated  angle-rate  is  then  fed  back  through  an  integrator  and  subtracted  from 
the  original  tracking  error  [19]. 

As  seen  from  Figure  2-5,  the  magnitude  of  the  nulling  angle-rate  signal  is 
proportional  to  the  tracking  error  represented  by  the  sloped  gain  curve.  Two 
curves  are  presented  in  Figure  2-5.  The  top  curve  represents  a  "loose"  tracking 
loop,  while  the  bottom  curve  represents  a  "tight"  tracking  loop.  Both  curves  show 
a  saturation  gain  which  represents  the  maximum  angle  rate  the  amplifier  can 


Figure  2-4b  Missile  Seeker  Tracking  Loop 


generate.  Below  the  saturation  gain  is  a  corresponding  angle  which  represents 
the  tracking  error  that  generates  saturation  (and  hence  the  maximum  tracking 
error  allowed  before  break- lock).  The  angle  resolution  of  the  amplifier  is  assumed 
to  be  between  5%  and  10%  of  the  saturation  angle.  The  amplifier  is  assumed  to  be 
updated  with  new  tracking  errors  at  the  seeker  scan-rate  of  100  Hz  [191. 

The  minimum  resolvable  tracking  angle  for  Figure  2-5a  is  taken  as  10%  of 
®error  max  or  0.1°.  The  minimum  resolvable  tracking  angle  for  Figure  2-5b  0error 
min,  is  taken  as  5%  or  0.005°.  The  two  tracking  errors  are  multiplied  by  their 

corresponding  amplifier  gains  to  generate  angle-rate  errors  of 

maximum  error-rate  (600/sec)  (1/200°)  =3°/sec  (2-26a) 

minimum  error-rate  (15/sec)  (1/10°)  =  1.5°/sec  (2- 26b) 

which  gets  passed  to  the  missile  autopilot  at  100  Hz.  Thus  a  good  range  for  a  Q(t) 
value  is  between  900°^/ sec^  and  225°^/ sec^  where  Q(t)  is  defined  as  the  white  noise 
diffusion  strength  [19]. 

The  second  noise  source  represents  the  missile  autopilot's  ability  to  produce 
identical  missile  accelerations,  given  identical  acceleration  commands  from  the 
seeker.  Figure  2-6  shows  a  reproduced  oscilloscope  plot  of  autopilot  acceleration 
error  versus  frequency  for  an  actual  missile  analyzed  at  the  Electronic  Warfare 
Division's  DIME  facility,  Wright  Laboratory,  Avionics  Directorate  [19].  The  plot 
shows  the  autopilot  transfer  function  to  be  an  approximate  first-order  lag  with  a 
300  Hz  cutoff  frequency  and  DC  gain  of  about  0.08  g's.  The  error  plot  can  be 
modelled  approximately  as 

autopilot  acceleration  error  0.08(g)  2000  /  (s  +  2000)  (2-27) 

Since  the  autopilot  error  bandwidth  is  much  greater  than  the  autopilot  transfer 
function  bandwidth  in  the  Cusumano  and  DePonte  model  (Figure  2-2),  the  autopilot 
noise  can  be  modelled  as  white  noise  of  strength  (0.08  g)2/second. 
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Figure  2.5a  Seeker  Amplifier  Gain  Curve  1.0  deg  Tracking  Error 


Figure  2.5b  Seeker  Amplifier  Gain  Curve  0.1  deg  Tracking  Error 
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Figure  2-6  Ajitopiloti  Acceleration  Ejror  Vs  Frequency 


2.4.2  Filter  Tuning  Noise 

The  filter  is  given  pseudo-noise  for  the  proportional  navigation  constant 
(n),  the  missile's  Surface  to  Mass  Ratio  (S/M),  and  the  estimated  scenario  altitude 
(h)  as  shown  in  Equations  (2-28  a-c) 


an/at  =  0.0  + wn(t) 

(2-28a) 

ah/at  =  0.0  +  w^t) 

(2-28b) 

a(S/M)/at  =  0.0  +  w(s/m)(*) 

(2-28c) 

The  pseudo-noise  is  added  to  compensate  for  initial  filter  uncertainties  for  the  two 
missile  parameters,  and  the  scenario  altitude.  Noise  is  also  added  to  the  filter’s 
lateral  acceleration  estimate  and  the  aircraft  acceleration  states.  The  missile  noise 
allows  the  filter  to  track  the  fifth-order  truth  missile.  The  aircraft  acceleration 
noise  allows  the  filter  to  track  the  pilot-induced  maneuver  errors  discussed  in 
Section  2.4.1.  All  noise  strengths  are  left  as  tuning  parameters  for  the  filter.  In 
tuning  the  filter,  the  noise  strengths  are  kept  as  small  as  possible  to  allow 
acquisition  and  track.  The  "low  noise"  tuning  philosophy  sacrifices  filter  settling 
time  to  achieve  a  better  final  volume  estimate  through  slower  diffusion  of  tuning 
noise  into  the  final  volume  estimate. 

2.4.3  Ladar  Measurement  Noise 

Ladar  measurement  noise  originates  form  three  sources:  Turbulence  noise 
affects  the  ladar's  detector  array  resolution  of  angle  and  angle-rate,  0(t)  and 
a0(t)/at,  measurement;  receiver  processing  delay  noise  which  degrades  range, 
R(t),  measurement;  and  spectral  frequency  resolution  which  degrades  range-rate 
measurement,  aR/at .  Noise  statistics  are  provided  in  Table  2-3. 

Figure  2-7  shows  the  structure  of  a  typical  ladar  tracking  system  system. 

The  primary  subsystems  include  the  laser,  the  receiver  optics,  the  infrared 
detector  array,  and  the  signal  processing  elements.  The  laser  transmits  a  pulse 
which  is  steered  to  the  expected  location  of  the  missile  through  use  of  beam  control 
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Table  2-3  Ladar  Measurement  Noise 


MEASUREMENT  TYPE 

RMS  Value 

Range 

0.5  meter 

Range-rate 

2  meter/sec 

Angle 

40.0  arad 

Angle-rate 

2.0  mrad/sec 

optics  and  slew  mirrors  (a).  The  laser  beam  divergence  is  set  to  equal  the  FOV  of 
the  detector  array.  The  missile  scattered  laser  return  is  focused  by  the  receiver 
optics  onto  the  detector  plane  (b).  When  the  received  pulse  crosses  a  set  detection 
threshold,  the  pulse  is  tagged  by  the  next  available  clock  pulse  from  the  receiver 
clock  (c).  Since  the  transmitted  pulse  is  released  by  a  signal  from  the  transmitter 
clock,  (presumed  to  be  in  synch  with  the  receiver  clock),  the  range  is  calculated 
from  the  elapsed  clock  time. 

Turbulence  effects  are  modeled  using  the  refractive  index  structure 
parameter  Cn2.  The  Cn2  parameter  provides  an  approximate  means  of  estimating  the 

rms  angle  error  caused  by  propagation  through  a  turbulent  atmosphere,  and  is 
given  by  Equation  (2-29)  [161. 

C^2  =  4.2  1C 14  h' 1/3  e-M®  (2-29) 

<02>=Cn29tDO'1/3  (2-30) 


where, 

C^2  is  the  refractive  index  structure  parameter, 

h  is  the  scenario  altitude, 
ho  is  3200  meters, 

91  is  the  target  range, 

<02>  is  the  mean  squared  beam  wander  angle, 
Dq  is  the  optic  diameter. 
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Figure  2-7  A  Typical  Ladar  System 


The  received  signal  is  picked  up  by  the  detector  array  (shown  in  Figure  2-8) 
which  maps  the  system  FOV  into  the  individual  detector  element  FOV’s.  An  angle 
error  signal  is  then  generated  by  the  receiver  which  enables  the  ladar  to  measure 
angle  and  angle  rate.  The  FOV  of  the  ladar  is  determined  by  the  number  of 
detectors  in  the  detector  array,  the  scan-rate  (the  ratio  of  the  ladar’s  FOV  to  the 
detector  array  FOV)  -  assumed  to  be  one  -  and  the  instantaneous  FOV  of  each 
detector  [Seyrafi]. 

Ladar  FOV  <xl1/2  =  n<j1/2  kscan1/2  na^2  (2-31) 

where 

«L1/2  is  the  FOV  plane  angle, 

n<j  is  the  number  of  detectors  (256), 

kscan  is  the  scan-rate  (assumed  to  be  one), 

ad1/2  is  the  detector  FOV  plane  angle  (40  urad  assumed). 

Normally  in  optics,  FOV  angles  are  given  as  solid  angles.  However  since  this  thesis 
is  limited  to  a  planar  analysis,  the  notation  a 1/2  is  used  to  represent  planar  angles. 
Since  the  FOV  of  the  ladar  is  small  (640  urad  under  the  assumed  conditions),  the 
ladar  depends  on  the  IR  staring  array  of  the  MWR  to  perform  low  resolution 
angular  tracking.  The  ladar  is  then  cued  in  to  allow  high  resolution  angle 
measurements. 

Figure  2-8  shows  part  of  a  16  x  16  element  detector  array.  Each  detector  maps 
into  one  part  of  the  receiver  FOV.  The  inner  solid  circle  represente  the  focused 
return  signal.  The  outer  ring  represents  the  one  sigma  error  bound,  <e2>l/2,  which 
determines  the  angular  resolution  of  the  ladar.  The  value  of  <02>^/2  has  been 
demonstrated  to  about  40  urad  in  tested  systems,  and  is  verifiable  using  Equations 
(2-29)  and  (2-30).  Because  a^l/2  of  each  detector  is  40  urad,  the  outer  ring  would 

normally  be  approximately  equal  to  the  individual  detector  size,  but  has  been 
enlarged  in  Figure  2-8  for  illustration  purposes. 
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Most  detector  noise  sources,  such  as  bandwidth  noise,  dark  current,  shot 
noise,  etc.,  affect  the  signal-to-noise  ratio  of  the  detector  and  hence  the  maximum 
effective  range  of  return  pulse  detection  [2].  For  a  typical  1.064  nm  Nd:YAG  pulsed 
laser,  the  maximum  detection  range  is  well  over  30  km  in  clear  atmosphere.  Since 
the  scenarios  studied  in  this  analysis  are  well  under  that  range  (typically  10  km  or 
less)  receiver  detection  noise  is  not  considered  applicable.  The  detector  noises  do 
not  affect  the  range  and  angle  resolution  of  the  receiver. 
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Most  receivers  use  a  silicon  based  detector  head  with  variable  but 
measurable  rise  times  of  between  2-50  ns.  The  rise  times  are  accurately  measured 
and  can  be  compensated  for  by  the  detection  software.  Therefore  detector  rise  time 
does  not  significantly  affect  range  resolution. 

When  a  laser  pulse  is  transmitted  and  received,  an  internal  receiver  clock 
keeps  track  of  the  transmit  and  receive  times  through  an  integer  count  of  the 
elapsed  clock  periods.  Therefore  the  Time  of  Launch  (Tqj  )  and  Time  of  Arrival 
(Tqa)  resolution  of  the  ladar  cannot  exceed  the  resolution  of  the  receiver  clock. 
Current  laser  trackers  can  support  a  Tqa  resolution  of  about  2.5  ns,  as  illustrated  in 
Figure  2-9  [211.  With  a  similar  Tqi  resolution,  this  gives  a  receiver  clock  induced 
range  error  of: 

range  error  9te  =  ±Tclocknair c  (2-32) 

where 

Tciock  is  the  receiver  clock's  resolution  (2.5  x  10'^  sec), 
nalr  is  the  index  of  refraction  of  air  from  Equation  (2-33), 

c  is  the  speed  of  light  in  a  vacuum  (2.998  x  10®  m/sec). 

From  Equation  (2-32)  the  rms  range  error  is  approximately  0.75  meters. 

Atmospheric  turbulence  is  another  potential  source  of  ladar  range 
measurement  noise.  The  speed  of  light  in  any  non- vacuum  medium  is  given  by  the 
speed  of  light  in  vacuum,  c,  divided  by  the  index  of  refraction  in  that  medium,  with 
the  index  of  a  vacuum,  n0,  defined  as  1.0.  According  to  Goodman,  the  index  of 
refraction  of  air,  nair,  varies  as  a  function  of  laser  wavelength,  X,  atmospheric 
pressure,  P,  and  absolute  temperature,  T  [7]: 

index  of  air  najr=  1  +  77.6(1  +  7.52  x  10"®  X'^)  P  x  lO’^/T  (2-33) 

where 

X  is  the  wavelength  in  jim  (1.064), 

T  is  the  atmospheric  temperature  in  K°> 

P  is  the  atmospheric  pressure  in  millibars. 
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Figure  2-9  Receiver  Timing  Diagram 


Turbulence  is  caused  by  small  scale  variations  in  atmospheric  temperature  and 
pressure.  The  net  effect  of  turbulence  is  to  cause  a  corresponding  variation  of  the 
index  of  refraction,  and  therefore  of  the  speed  of  propagation.  Because  range  is  a 
large  scale  measurement,  the  effective  index  of  refraction  will  be  an  average  of 
the  total  of  all  refraction  indexes  traversed  appropriately  scaled  to  account  for 
propagation  distances  crossed  with  that  index.  When  the  average  effect  is 
accounted  for,  turbulence  has  a  negligible  effect  on  range  error  at  the  ranges  of 
interest.  This  argument  is  made  on  three  premises:  That  index  variation  is  a  linear 
function  of  pressure,  that  the  index  only  varies  slightly  with  pressure,  and  that  the 
total  pressure  variations  encountered  must  average  to  approximately  the  nominal 
atmospheric  pressure.  When  all  three  premises  hold,  the  net  range  error  caused 
by  atmospheric  turbulence  is  less  than  10'^  x  total  distance  propagated,  or  less  than 
0.1  meter  per  10  km  of  tracking  range.  Since  the  scenarios  encountered  in  this 
study  are  usually  simulated  at  ranges  under  10  km,  the  effect  of  turbulence  on 
range  measurement  is  negligible. 

The  ladar  range-rate  noise  can  be  calculated  using  Equation  (2-34)  f  16:237] 

Af  =  2  (d9t(t)/dt)  /  k  (2-34) 

where 

Af  is  the  Doppler  frequency  shift, 
d9t(t)/at)  is  the  range-rate, 
k  is  the  laser  wavelength, 

by  setting  the  frequency  shift  equal  to  the  frequency  resolution  of  a  heterodyne 
receiver.  (A  heterodyne  receiver  uses  a  frequency  mixer  to  generate  a  beat 
pattern  between  a  received  signal  and  a  reference  signal.  The  frequency  of  the 
beat  pattern  determines  the  frequency  difference  between  the  received  signal  and 
the  reference  signal.)  A  1.0  m/sec  range-rate  resolution  requires  a  receiver 
frequency  resolution  of  approximately  2.0  MHz.  Seyrafi  verifies  that  range-rate 
accuracies  less  than  1.0  m/sec  are  achievable  [16:238]. 
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2.5  The  Initial  Conditions 


Initial  conditions  are  established  by  setting  up  a  particular  scenario  of 
interest.  Initial  missile  range  and  tracker  angles  are  chosen,  along  with  missile 
velocity,  aircraft  heading,  aircraft  turning  rate,  the  proportional  navigation 
constant,  the  missile's  surface-to-mass  ratio,  and  scenario  altitude.  The  specific 
initial  conditions  chosen  for  these  variables  are  left  to  the  problem  parameter 
under  investigation. 

The  initial  conditions  are  assumed  to  be  handed  off  from  a  Missile  Warning 
Receiver  (MWR)  onboard  the  aircraft.  Initial  conditions  to  be  provided  include 
range,  range-rate,  time-to-go,  and  LOS  angle.  Values  for  the  rms  errors  of  the 
initial  conditions  are  based  on  minimum  performance  goals  established  by  Mr 
Kevin  McCamey,  Program  Manager  of  the  Electronic  Countermeasures  Advanced 
Development  Branch  Missile  Warning  Program!  12],  and  are  provided  in  Table  (2-4). 

Initial  covariance  estimates  for  the  aircraft  states  are  assumed  to  be 
available  from  the  GPS/INS.  GPS/INS  covariance  values  are  obtained  from  Negast 
[13]  using  the  twenty-two  state  GPS/INS  model  in  steady  state.  Table  (  2-5  )  contains 


Table  2-4  MWR  Provided  Initial  Conditions 


Initial  Parameter 

RMS  Error 

Range 

3.0  meters 

Range  Rate 

5.0  meter/ sec 

Angle 

2.0mrad 

Time-to-go 

10%  error 
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the  rms  errors  associated  with  the  aircraft  states,  rhese  same  rms  errors  are  used 
in  the  measurement  update  cycles  tor  the  aircraft  states,  since  they  are  assumed  to 
be  continuously  available  from  the  navigation  system. 


fable  2-5  INS  Provided  Initial  Conditions 


INS  State 

RMS  Error 

North  Velocity 

0.03  meter/sec 

Fast  Velocity 

0.03  meter/ sec 

Vertical  Position 

9  feel 

The  filter's  proportional  navigation  parameter,  n,  is  given  an  initial 
covariance  of  0.25  guaranteeing  that,  with  an  initial  estimate  of  4.5,  the 
parameter's  full  range  (  3  to  6  )  is  covered  within  the  3o  bounds.  The  filter's  missile 
surface-to-mass  ratio  is  given  a  covariance  of  (0.0008)2  to  cover  the  expected  range 
(0.005  to  .0001)  with  an  initial  estimate  of  0.0025  m2/kg. 

2.6  The  Summarized  Truth  and  Filter  Model  Equations 

This  section  presents  the  actual  model  equations  and  noise  matrices  used  in 
simulation.  A  rigorous  presentation  of  extended  Kalman  filtering  is  not  presented, 
nor  are  the  equations  associated  with  extended  Kalman  filtering.  Knowledge  of  the 
subject  is  assumed  on  the  part  of  the  reader.  For  a  rigorous  treatement  of  extended 
Kalman  filtering  refer  to  May  beck  [9]. 

The  system  dynamics  are  developed  in  Section  2.2,  the  state  selections 
summarized  in  Tabie  2-2,  and  the  process  and  measurement  noise  terms  presented 
in  Section  2.4.  The  following  is  a  summarized  list  of  the  actual  truth  model  and 
filter  model  equations  used  in  simulation.  All  variables  and  notation  used  here  is 
consistent  with  that  used  throughout  the  chapter.  The  states  are  as  defined  in 
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Table  2.2a  and  2-2b.  The  truth  and  sytem  models  are  described  by  two  general 
equations  f9|: 


ax/  at  =  ifx(t),t)  +  G(t)w(t)  (2-35) 

z(t,)  =  lifx(ti),til  +  y(tj)  (2-36) 

where 

f[x(t),t)  is  a  vector  of  non-linear  state  equations  describing  system  dynamics, 
>£(t)  is  a  vector  of  states  from  either  Table  2-2A  or  Table  2-2B, 
w(t)  is  a  vector  of  dynamic  process  noises, 

G(t)  is  a  non-linear  noise  shaping  matrix, 
z(tj)  is  a  vector  of  system  measurements, 
il[x(t|),tj|  is  a  vector  of  non-linear  measurement  equations, 
v(tj)  is  a  vector  of  measurement  noises  or  uncertainties. 

2.6.1  The  Truth  Model  Dynamics 

Equation  bets  (2-37)  and  (2-38)  define  the  true  dynamics  vector  a^/at: 


3Xsi/3t  =  -nvm1/2xsl  -  f2Aj  2  4-  A)  xs2/vm  (2~37a) 

0xs2/dt  =  -  flVm1/2xs2  -  f2AL2  -  Al  Xsl/Vm  (2-37b) 

3xS3/3t  =  xs4  (2-37c) 

3x54/ at  =  \  -  AtSinfxS3  -  ArcTan(x5]Q/xS9)]  +  (2-37d) 

AmSin(xS3  -  ArcTan(3X52/3t  /  axjj/at)]  +X56Xs4  J  /\^5 

dxs5/dt  =  xs6  (2-3  7e) 

dXs5/dt  =  AtCos[Xs3  -  ArcTanlXjjo/Xsg)]  -  (2-37f) 

AmCosfXs3  -  ArcTan(axs2/3t  /  ax^/atM  +  Xs5Xs42 

dX57/dt=X59  (2~37g) 

axs8/at  =  xsl0  (2-37h) 

3Xs9/at=  to^RturnSin(o)t+e)  +  Wg(t)  (2-37i) 

axsio/at  =  -  <03RturnCos(«Jt+e)  +  W10(t)  (2-37j) 

dXji  i/3t=  -  13.33x5! !  +  13.33x54  +  W|  j(t)  (2-37k) 

dX5]2/at  =  -3.33xsi2  +  3.33nx5j  jxs6  (2-371) 

dXsi3/dt=  -  3..06X5J3*  4( .99875  1)x512  +  w13(t)  (2-37m) 


<1*5 14/ at  =  -  1 4  -  4(.783663)xsi  2  (2-37n) 

axsl5//at=  -  3..0bx<5  ]  5  -  4(  1 .1 208(>)xs  j  2  (2-37o) 

and 

Vm=(xsi2  +  xs22)1/2  (2-38a) 

Am  =  (axsl/at2  +<)xS2/(n2)l/2  (2-38b) 

At  =  (Xs92  +  x^)2) 1/2  (2 -38c) 

Ai,  =  xsl3  +  xsl4  +  xsl5  (2-38d) 

f,  =  1.29e'(h/30>000)  S/My  1/2  (2-38e) 

f2  =  2.0/(1.29  e'<h/30-000>  S/My)  (2-380 


2.6.2  The  Truth  Model  Measurements 

Equation  Set  (2-39)  defines  the  true  measurement  vector  ^(t,): 

Zsl  =  -  (Xs6+vs6)tos(Xs3+vs3)  +  (Xs5+vs5HXs4+vs4)sin<xs3+vs3)  (2-39a) 
+  <Xs7+vs7) 

Zs2  =  -  (Xs6+vs6>Sin<xs3+vs3>  '  (xs5+vs5)(xs4+vs4)Cos(xs3+vs3)  (2-39b) 


Zs3  =  *s3  +  vs3  (2-39c) 

2^4  =  Xs4  +  vs4  (2-39d) 

Zs5  =  xs5  +  Vs5  (2-39e) 

Zs6  =  xs6  +  Vs6  (2-39f) 

zS7  =  xS7  +  vS7  (2-39g) 

zs8  =  xs8  +  vs8  (2-39h) 

Zs9  =  Xs9  +  vs9  (2-39i) 

zslO==xslO  + vslO  (2-39|) 


2.6.3  The  Truth  Model  Noise  Strength 

Table  2-6  presents  the  noise  strength  values  used  in  the  dynamic’s  Q(t) 
matrix  and  in  the  measurement's  R(tj )  matrix.  Each  noise  strength  is  presented 


2-33 


with  its  accompanying  subscripted  rm  or  qnn  term  to  identify  matix  location.  All 
matrix  entries  not  explicitly  specified  in  the  table  are  assumed  to  be  zero. 


f  able  2-6  Truth  Model  Noise  Terms 


R  Entry 

Strength 

Q.  Entry 

Strength 

r3 

1.6  E" 10 

49.9 

1.0 

r 4 

4.0  F  6 

410.10 

1.0 

rS 

0..25 

4ll. ll 

15.0 

1-6 

1.0 

413.14 

.64 

r7 

.01 

r8 

0.01 

r9 

25.0  F 10 

rio 

25.0  F 10 

2.6.4  The  Filter  Model  Dynamics 

Equation  Sets  (2-40)  and  (2-41)  define  the  filter  dynamics  vector  ax^/at: 


axfl/dt  =  -  flVm1/2xn  -  f2AL2  +  ALxf2/Vm  (2-40a) 

axf2/at  =  -  flVm1/2xf2  -  f2AL2  -  ALxfl/Vm  (2-40b) 

axf3/at  =  xf4  (2-40c) 

axf4/at  =  { -  At  Sin[xf3  -  ArcTan(xfl5  /  xn3)]  +  (2-40d) 

AmSin[xf3  -  ArcTan(axf2/at/  ax^/at)]  +  xffcXf4  }  /  xf5 

axf5/at=xf6  (2-40e) 

axf6/at  =  At  Cosfxf3  -  ArcTan(xf15  /  xfl3)J  -  (2-40f) 

Am  Cosfxf3  -  ArcTan(axf2/at  /  axf^at))  +  xf5xf42 

axfy/at  =  0.0  (xj7  =  1 .5xf5/xf6)  ( 2-40g) 

axf8/at  =  0.0  + wf8(t)  (2-40h) 

axf9/at  =  0.0  +  wfg(t)  (2-40i) 
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3Xfio/at  =  0.0  + wno(t)  (2-40j) 

dxn  \/dl  =  0.0  + wt  ]  i(t)  (2-40k) 

«9xf 1 2/at  =  xn  3  (2-401) 

dXf ,  3/dt  =  w  3R  lurnSi  n  («>t+«)  +  wfl3(t)  (2-40m) 

OXf}4/dl=  xf15  (2-40n) 

axfl5/at  =  ‘  «'>3RturnCos(wl+fl)  +  wfl5(t)  (2-40o) 

and 

vm  =  (xfl2  +xf22)1/2  (2-41a) 

Am  =  (axfi/at2  +  oxQ/at2) 1/2  (2-4  lb) 

At  =  (xfl32  +  xfl52>1/2  (2-41c) 

Al  =  n  Xf4  Xf6  +  xf ,  o  (2-41d) 

f !  =  1.29  e'(xfl  1/30.000)  Xfg  y  1/2  (2-41e) 

f2  =2.0/(1.29  e’<xn  1/30.000)  Xfgy)  (2-410 


2.6.5  The  Filter  Measurement  Matrix 

Equation  Set  (2-42)  presents  the  filter  measurement  vector  z^tj): 


zfl  =  '  (xf6+vf6)Cos<xf3+vf3)  +  <xf5+vf5 > < xf4+vf4) Sin ( xf3+vf3 )  (2-42a) 

+  (xfl2+Vf12) 

zf2  =- (xf6+vf6)Sin(xf3+vf3)  -  (Xf5+vf5)(xf4+Vf4)Cos(xf3+vf3)  (2-42b) 
+  (Xf14+vfl4) 

zf3  =  xf3  +  vf3  (2-42c) 

Zf4  =  Xf4  +  vf4  (2-42d) 

zf5  =  xf5  +  vf5  (2-42e) 

zf6  =  xf6  +  vf6  (2-42f) 

Zfi2  =  Xfi2  + Vfi2  (2-42g) 

zfl3  =  xfl3  +  vfl3  (2-42h) 

zfl4  =  xfl4  +  vfl4  (2-42i) 

zfl5  =  xfl5  +  vfl5  (2-42j) 
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The  filters  measurement  matrix  (H)  is  formed  by  taking  the  partial 
derivatives  of  Equations  (2-42  a-j)  with  respect  to  the  15  filter  states.  The  11  matrix 
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where 


a  =  Xf6Sin(Xf3)  +  Xf5Xf4Cos(xf3) 

(2-43a) 

b  =  Xf5Sin(Xf3) 

(2-43b) 

c  =  Xf4Sin(Xf3) 

(2-43c) 

d  =  Cos(Xf3) 

(2-43d) 

e  =  -  Xf5Cos(xf3>  +  Xf5Xf4Sin(Xf3) 

(2-43e) 

f  =  -  xf5<J°s(xf3) 

(2-430 

g  =  -  Xf4Cos(Xf3) 

(2-43g) 

h=  -  Sin(xf3) 

(2-43h) 

2.6.6  The  Filter  Model  Noise  Strength 

Table  2-7  presents  the  noise  strength  values  used  in  the  filter  dynamic’s  Q;  *4 
matrix  and  in  the  filter  measurement's  R(tj)  matrix.  Each  noise  strength  is 
presented  with  its  accompanying  subscripted  r^  or  term  to  identify  matix 
location.  All  matrix  entries  not  explicitly  specified  in  the  table  are  assumed  to  be 
zero. 
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l  able  2-7  Filter  Model  Noise  Terms 


R  Entry 

Strength 

Q.  Entry 

Strength 

rl.l 

aU 

48.8 

.25 

r2.2 

b2.2 

49.9 

1.0 

r3.3 

1.6  r 10 

410.10 

64.0 

r4.4 

4.0  f  6 

4ll.ll 

1.0  E- 8 

r5.5 . 

025 

413.13 

1.0 

r6.6 

1.0 

415.15 

1.0 

r7.12 

.01 

r8.l3 

25.0  F7 10 

r9.14 

.01 

r  10.15 

25.0  F7 10 

where 

al,l  =  r4i4Xf52Sin(xf3)  +  r6)6Cos(xf3)  (2-44a) 

b2,2  =  r4(4xf52Cos(xf3)  +  r6>6Sin(  xf3)  (2-44b) 


2.7  Summary 

Chapter  II  developed  the  truth  model  and  filter  model  used  in  the  simulation 
studies  conducted  for  this  thesis.  Section  2.1  defined  the  simulation  reference 
irames  and  gave  rationale  for  their  selection.  Section  2.2  developed  the  kinematic 
equations  describing  aircraft  and  missile  motion.  Section  2.3  described  the 
selection  of  truth  model  and  filter  states  based  on  the  equations  developed  in 
Section  2.2.  Section  2.4  presented  the  noise  sources  driving  the  aircraft  dynamics, 
missile  dynamics,  and  ladar  measurements,  and  gave  rationale  for  the  noise 
strengths  used  in  the  simulation  models.  Section  2.5  covered  the  selection  of  initial 
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conditions  which  were  considered  to  have  special  relevance  to  the  simulation 
model.  Section  2.6  summarized  the  state  equations,  measurement  equations  and 
noise  matrices  used  in  the  simulation.  Hie  next  chapter  overviews  the  simulation 
runs  conducted  in  this  study,  bach  run  focuses  on  an  individual  parameter  viewed 
to  be  significant  in  determining  filter  either  performance  or  truth  model 
behavior. 
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Ill  The  Simulation  Runs 


3.0  Introduction 

This  chapter  presents  an  overview  of  the  computer  simulations  conducted 
while  researching  the  kinetic  kill  zone  estimation  problem.  All  computer 
simulations  presented  in  this  thesis  study  were  accomplished  on  SUN  SPARC  II 
workstations,  using  the  MSOFE  simulation  software  (3],  and  standard  FORTRAN  77 
computer  code.  Chapter  III  only  presents  parameters  varied  in  each  run  (and  the 
values  used  in  simulation).  Actual  simulation  results  are  presented  in  Chapter  IV. 
Section  3.1  presents  the  missile-aircraft  geometry  run  series.  Eight  runs  were 
made  to  analyze  the  effect  of  the  missile  approach  angle  on  the  ability  to  predict 
the  kill  zone.  Section  3.2  presents  the  ladar  run  series,  conducted  to  analyze  how 
ladar  measurement  parameters  affect  the  ability  to  predict  the  kill  zone.  Each 
measurement  is  first  varied  independently,  then  all  measurements  are  varied 
simultaneously,  and  finally  the  ladar  measurement  rate  is  varied  from  60  Hz  to  30 
Hz.  Section  3.3  presents  a  run  series  conducted  to  analyze  how  specific  aircraft 
parameters  affect  the  final  predicted  volume.  Section  3.4  presents  a  worst  case 
performance  baseline.  (The  best  case  performance  baseline  is  taken  from  Section 
3.1.) 

As  a  final  note,  certain  terminology  consistently  used  throughout  this 
Chapter  needs  a  brief  clarification.  The  term  "run  series"  always  refers  to  a 
collective  set  of  all  runs  done  to  analyze  one  of  three  major  subgroups  affecting 
kill  zone  prediction:  missile-aircraft  geometry,  the  ladar,  or  aircraft  specific 
parameters.  The  term  "run  set'  always  refers  to  a  sequence  of  runs  made  to  analyze 
one  specific  parameter  within  a  run  series.  The  term  "run,"  by  itself,  always  refers 
to  a  single  computer  simulation. 
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3.1  The  Missile  Approach  Angle  Run  Series 

The  missile  approach  angle  simulations  presented  here  are  conducted  to 
analyze  the  impacts  of  missile  engagement  geometry  on  the  final  volume 
predictions  achieved.  For  all  missile  approach  angle  simulations,  the  aircraft  is 
assumed  to  be  initially  heading  east  at  300  m/s  and  starting  a  3  g  turn  to  the  north. 
(The  aircraft  maneuver  chosen  as  a  baseline  condition  against  which  a 
comparative  analysis  is  later  performed.  The  airspeed  and  turn-rate  were  selected 
as  being  a  typical  mid-range  value.)  Missile  engagements  are  taken  in  ji/4 
intervals  starting  at  0.0  radians  and  measured  in  the  missile-to-aircraft  reference 
frame  (i.e.,  the  initial  LOS  is  defined  looking  from  the  missile  to  the  aircraft). 

Initial  missile  speed  is  1500  m/s  for  all  aspect  angles.  Table  3-1  summarizes  the 
initial  conditions  used  in  the  eight  missile  approach  angle  runs. 

Due  to  the  extensive  variations  in  missile  kinematics  which  results  from  the 
different  aspect  angles,  there  are  significant  differences  in  two  missile  initial 
conditions  other  than  approach  angle.  The  two  variable  initial  conditions  are 
range  and  missile  surface-to-mass  ratio.  The  variations  in  these  missile  initial 
conditions  are  driven  primarily  by  the  dynamics  of  a  post-bumout  coasting 
missile.  In  each  aspect  angle,  the  set  of  initial  conditions  which  resulted  in  the 
best  missile  kinematics  ("best"  from  a  simulation  viewpoint,  not  necessarily  from  a 
cross-section  prediction  viewpoint)  is  ultimately  chosen. 

The  missile's  surface  to  mass  ratio,  S/M,  is  set  near  the  minimum  of  the 
normal  expected  range  (0.0001  to  0.005  m2/kg)  for  the  n/4  approach  angle  because 
ic/4  is  a  very  benign  dynamic  approach  angle,  and  the  missile  is  initially  chasing 
the  aircraft.  Reducing  S/M  increases  the  coasting  performance  of  the  missile  and 
therefore  allows  a  longer  missile  flight  time.  The  increased  range  and  missile 
flight  time  allows  the  final  kill  zone  volume  to  be  more  easily  analyzed  as  a 
function  of  range  and  simulation  time  (see  Chapter  IV).  As  the  initial  missile  Table 
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4-1  Missile  Approach  Angle  Series 


approach  changes  from  tail-chase  to  head-on  approach,  the  missile's  S/M  also 
increases  to  allow  the  missile  better  turning  performance  at  the  expense  of  poorer 
coasting  performance.  (Refer  to  Chapter  II  Subsection  2.2.4  for  the  relationships 
between  S/M  and  missile  drag).  All  missile  approach  angle  trajectory  plots,  filter 
plots  and  cross-section  prediction  plots  are  presented  in  Appendix  A. 

There  are  two  aspect  angles,  3n/4  and  7n/4,  in  which  the  simulation  fails 
altogether.  ).  The  simulation  results  indicate  that  the  missile  is  unable  to  maintain 
track  on  the  aircraft  for  the  3n/4  and  7n/4  approach  angle's,  so  these  geometries 
are  presented  now  and  are  not  included  with  the  results  presented  in  Chapter  IV. 

Figure  3-1  shows  a  vector  diagram  explaining  why  the  3*/4  (and  by  analogy 
the  7re/4)  missile  shot  dynamics  failed.  The  missile  commands  lateral  acceleration, 
as  given  in  Equation  (2-  13),  proportional  to  the  LOS-rate  in  the  missile-to-aircraft 
reference  frame.  As  seen  from  Figure  3-1,  there  are  three  factors  which 
contribute  to  the  total  LOS-rate:  drag  induced  LOS-rate,  range-closure  induced 


3-3 


trajectory  after 
3g  turn  begins 

d 


aircraft  initial 
trajectory 


missile  north 
velocity 


A 


KEY 


A  is  the  initial  missile-to-aircraft  LOS 
Bl  is  the  initial  missile  trajectory  N. 

B2  is  the  missile  trajectory  after  drag  losses^ 
b  is  the  aircraft  initial  trajectory 
B 1  -  C  is  the  LOS  change  due  to  drag 
d  is  the  aircraft  trajectory  after  3g  turn 
D  -  Bl  is  the  LOS  change  due  to  aircraft  turn 
E  is  the  final  missile-to-aircraft  LOS 
E  -  D  is  the  LOS  change  due  to  range  closure 
F  is  the  total  LOS  change  due  to  all  effects 


missile  east 
velocity 


Missile  Seeker  (5iigin 


Figure  3-1  Failed  Missile  Kinematics 
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LOS- rate,  and  aircraft  turn  induced  LOS-rate.  The  scenario  starts  with  an  aircraft 
trajectory  defined  by  b,  a  missile  trajectory  defined  by  B 1,  and  an  LOS  defined  by 
A.  If  no  aircraft  and  missile  trajectory  changes  occur,  the  LOS  vector,  A,  never 
changes  and  the  missile  intercepts  the  aircraft  When  missile  drag  is  added  to  the 
scenario,  the  missile's  trajectory  changes  to  B2,  and  a  drag  induced  LOS-rate  occurs 
given  by  the  angle  B  \  -  C.  When  the  aircraft  turn  is  added  to  the  scenario,  the 

aircraft  trajectory  changes  from  b  to  d  and  the  turn  induced  LOS-rate  is  given  by 
the  angle  D  -  Bj.  When  range  closure  is  considered  a  third  source  of  LOS-rate  is 

given  by  the  angle  E  -  D.  The  combined  angle-rate  is  given  by  the  arc  F.  In  the 
3n/ 4  scenario,  the  drag  induced  LOS-rate  is  acting  against  the  turn  induced  LOS 
rate.  When  the  missile  commands  lift  to  null  the  aircraft's  turn  induced  LOS  rate, 
the  missile's  lift-drag  acts  to  increase  the  missile-to-aircraft  LOS-rate.  The  harder 
the  missile  turns,  the  harder  the  missile  needs  to  turn.  This  effect  is  heightened  by 
the  range  closure,  which  is  severe  for  a  forward  quarter  missile  approach.  The 
missile  literally  turns  itself  to  death. 

3.2  The  Ladar  Noise  Run  Series 

The  ladar  simulations  presented  here  were  conducted  to  analyze  the  impacts 
of  the  assumed  ladar  measurement  parametric  s  on  the  results  achieved.  An  ladar 
noise  run  series  are  simulated  from  the  n/4  missile  approach  angle,  shown  in 
Figure  3-2,  using  identical  initial  conditions  and  random  number  seeds  for  each 
run.  The  only  parameter  changed  for  each  run  is  the  specific  parameter  under 
investigation.  Table  3-2  shows  the  values  used  for  the  ladar  measurement  noise 
during  the  ladar  run  series.  The  baseline  values  are  selected  based  on  the  ladar's 
rms  errors  as  developed  in  Chapter  II.  The  other  values  are  generated  by 
redoubling  the  baseline  values  as  one  crosses  the  table  from  left  to  right.  The 
angle-rate  measurement  is  only  doubled  twice  since  it  is  believed  that  an  angle- 
rate  measurement  error  of  16.0  mrad/s  is  not  of  practical  value. 
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Table  3-2  RMS  Ladar  Noise 


Measurement 

Baseline  Value 

2nd  Value 

3rd  Value 

4th  value 

Range 

0.5  meters 

1.0  meters 

2.0  meters 

4.0  meters 

Range-rate 

1.0  m/s 

2.0  m/s 

4.0  m/s 

8.0  m/s 

Angle 

*40  urad 

80  ix  rad 

160  urad 

250  urad 

Angle-rate 

2.0  mrad/s 

2.0  mrad/s 

4.0  mrad/s 

8.0  mrad/s 

Table  3-3  presents  a  summary'  of  the  ladar  series  runs.  The  first  column  lists  each 
run  set  contained  in  the  ladar  series.  The  second  column  gives  the  parameter 
under  study  in  that  run-set.  The  values  used  under  the  "Parameter  Varied"  column 
are  listed  in  Table  3-2.  All  parameters  not  listed  in  the  "Parameter  Varied"  column 
are  left  at  the  baseline  values  (column  two.  Table  3-2)  for  the  entire  run  set.  The 
entry  "all"  indicates  that  all  parameters  listed  in  Table  3-2  are  varied 
simultaneously  from  the  baseline.  The  baseline  laser  pulse-rate  is  60  Hz,  with  a  run 
set  accomplished  at  30  Hz  for  comparison  as  requested  by  the  thesis  sponsor  [20]. 
The  last  column,  "Filter  vs.  System"  indicates  differences  between  the  filter’s 
assumed  rms  measurement  noise,  and  the  actual  rms  measurement  noise 
corrupting  the  system.  The  entry  "same"  indicates  that  the  truth  model's  and  filter 
model's  ladar  rms  measurement  noise  are  identical  for  the  parameter  under  study. 
The  entry  "filter  high"  indicates  that  the  filter's  noise  strength  is  varied  while  the 
system  noise  strength  remains  at  baseline.  The  entry  "filter  low"  indicates  that  the 
filter  noise  strength  is  kept  at  baseline  while  the  system  noise  is  varied. 

As  is  evident  from  Table  3-2  and  3-3,  the  ladar  series  addresses  five 
objectives:  (1)  Each  ladar  measurement  is  analyzed  independently  to  determine  the 
effect  of  rms  noise  on  the  final  volume  estimate;  (2)  All  ladar  measurements  are 
analyzed  as  a  group  to  determine  the  effect  of  ladar  measurement  quality  on  the 
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Table  3-3  Ladar  Series  Run  Sets 


Run  Set 

Parameter  Varied 

Filter  vs.  System 

Range 

Range  only 

same 

Range-rate 

i@Sii 

same 

Angle 

Angle  only 

same 

Angle-rate 

mmmm 

same 

60  Hz 

all  simultaneously 

same 

Range-f 

Range  only 

filter  tuned  high 

Range-s 

Range  only 

filter  tuned  low 

Range-rate-f 

filter  tuned  high 

Range- rate-s 

WKBBHBM 

filter  tuned  low 

Angle-f 

Angle  only 

filter  tuned  high 

Angle-s 

mmmm 

filter  tuned  low 

Angle-rate-f 

filter  tuned  high 

Angle- rate-s 

WKMBBBKSM 

filter  tuned  low 

all-f 

all  simultaneously 

filter  tuned  high 

all-s 

all  simultaneously 

filter  tuned  low 

30  Hz 

all  simultaneously 

same 

final  volume  estimate;  (3)  The  ladar  sampling  rate  is  analyzed  to  determine  its 
impact  on  the  final  volume  estimate;  (4)  The  ladar  is  tested  against  cases  where  the 
filter  assumed  measurements  are  more  accurate  than  the  true  measurements 
actually  are;  (5)  The  ladar  is  tested  against  cases  where  the  filter  assumed 
measurements  are  less  accurate  than  the  true  measurements  actually  are.  Since 
the  turbulence  model  used  throughout  this  thesis  study  assumes  benign 
turbulence,  the  effects  of  this  assumption  need  to  be  explored.  The  filter  "high”  and 
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filter  "low"  run  sets  are  conducted  to  determinetwo  things:  ( 1 )  Does  over  modelling 
turbulence  noise  in  the  filter  seriously  degrade  cross-section  estimates/  and  (2) 
Does  undermodeling  turbulence  noise  in  the  filter  cause  the  estimated  cross- 
section  to  not  capture  the  true  missile  trajectory  as  it  should?  All  ladar  run  series 
plots  are  presented  in  Appendices  B  (range),  (J  (range-rate),  D  (angle),  E  (angle- 
rate)  and  F  (pulse-rate). 

3.3  The  Aircraft  Run  Series 

The  aircraft  run  series  has  three  objectives:  ( 1 )  To  study  how  the  assumed 
DGPS/INS  accuracies  affect  the  final  kill  zone  volume;  (2)  To  study  how  the  aircraft 
turn  rate  affects  the  final  kill  zone  volume;  and  (3)  To  study  how  pilot-induced 
maneuver  error  affects  the  final  kill  zone  volume.  To  accomplish  these  objectives, 
one  run  set  is  generated  for  each  parameter  under  study,  using  the  n/4  geometry 
with  the  same  random  number  seeds  and  initial  conditions  as  in  the  ladar  series. 
The  simulation  model  originally  presented  in  Figure  1-4  assumes  that  acceleration 
measurements  can  be  provided  from  the  aircraft's  navigation  subsystem.  The 
DGPS/INS  run  set  removes  all  aircraft  acceleration  measurements  to  determine  the 
impacts  of  not  having  acceleration  measurements  available,  and  also  varies  the 
aircraft’s  velocity  measurement  noise  strengths  (as  shown  in  Table  3-4)  to 
determine  how  the  naivigation  subsystem's  estimated  errors  affect  the  predicted 
cross-section.  The  DGPS/INS  filter  plots  and  volume  plots  are  contained  in 
Appendix  H.  The  aircraft  turn  rate  set  is  generated  by  resetting  the  DGPS/INS  to 
the  normal  values  presented  in  Table  2-5,  and  varying  the  aircraft's  assumed  turn 
rate  from  0.0  g  to  5.0  g  in  1.0  g  increments.  The  aircraft  turn  rate  filter  plots  and 
volume  plots  are  contained  in  Appendix  I.  Finally  the  pilot  induced  maneuver 
error  is  analyzed  by  varying  the  aircraft  process  noise  strength  (Section  2-4)  from 
0.0  to  (0.5g)2/s  in  (0.1g)2/s  increments.  The  pilot  induced  maneuver  error  plots 
are  contained  in  Appendix  G. 
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Table  3-4  Assumed  DGPS/INS  Velocity  F.rrors 


Run 

North  Velocity  Error 

Oast  Velocity  Error 

1 

0.1  meier/second 

0.1  meter/second 

2 

0.5  meter/second 

0.5  meter/second 

3 

1.0  meter/ second 

1.0  meter/second 

4 

5.0  meter/second 

5.0  meter/second 

3.4  High  Noise  Baseline  Run 

One  final  run  is  made  to  establish  a  worst  case  volume  baseline.  The  "high 
noise"  baseline  is  established  by  taking  the  worst  assumed  noise  and  measurement 
conditions,  from  each  run  set,  and  using  those  conditions  together  in  a  single  run. 
The  ji/4  geometry  with  a  60  Hz  ladar  is  selected  for  consistency.  Table  3-5 
summarizes  the  lo  noise  strengths  used  to  establish  the  high  noise  baseline.  The 
plots  are  discussed  in  Chapter  IV. 


3.5  Summary 

This  Chapter  summarizes  the  MSOFE  simulation  runs  conducted  in  analyzing 
the  kinetic  kill  zone  prediction  problem.  The  parameters  analyzed  are  divided  into 


Table  3-5  High  Noise  Baseline  Run  Conditions 


Ladar 

Measurement 

lo  Noise  Strength 

Aircraft  Parameter 

lo  Noise  Strength 

Range 

4.0  meters 

North  Velocity 

5.0  m/s 

Range-rate 

8.0  m/s 

5.0  m/s 

Angle 

250  urad 

Turn-rate 

5.0  g 

Angle-rate 

8.0  mrad/s 

■HBH 

0.5  g 
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three  major  groups  -  or  series.  The  aircraft-missile  geometry  is  discussed  in 
Section  3.1.  The  iadar  measurement  parameters  are  covered  in  Section  3.2.  The 
aircraft  parameters  are  covered  in  Section  3.3.  The  x/4  missile  approach  angle 
geometry  is  used  as  a  best  case  performance  baseline  for  the  run  series'  conducted 
in  Section  3.2  and  Section  3.3.  A  worst  case  performance  baseline  is  developed  for 
the  same  run  series’  in  Section  3.4  using  the  worst  case  noise  strengths  from  each 
run  set. 

The  next  chapter  presents  the  simulation  results.  For  the  missile  appraoch 
angle  run  series  the  results  are  compared  to  analyze  actual  versus  expected  results. 
Two  run  sets  are  selected  from  the  missile  approach  angle  run  series  to  analyze  the 
effect  of  true  range  on  the  predicted  cross-section.  For  the  remaining  run  series, 
baseline  performance  criteria  and  volume  cross-section  figures  of  merit  are 
established  as  performance  comparison  tools. 


■ 
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IV  Results  and  Analysis 


4.0  Introduction 

This  chapter  presents  results  lor  the  simulation  series  discussed  in  Chapter 
111.  lhe  chapter  structure  parallels  the  structure  in  Chapter  111,  with  a  subsection 
designated  for  each  major  run  series.  Each  subsection  presents  the  results 
obtained  from  that  run  series,  a  brief  summary  of  the  results  and,  where 
applicable,  expected  results.  Discussion  focuses  in  three  primary  areas, 
explanation  of  the  results,  the  results'  effect  on  the  predicted  volume  cross-section, 
and,  where  applicable,  a  discussion  on  deviations  from  the  expected  results. 

Filter  and  volume  cross-section  plots  for  the  ji/4  missile  approach  angle  are 
presented  in  this  chapter  since  that  geometry  is  used  as  the  baseline  for  all  run 
series  except  the  missile  approach  angle  run  series.  Other  filter  plots  and  cross- 
section  estimation  plots  are  contained  in  the  appendices,  and  are  referenced  as 
needed. 

4.1  Missile  Approach  Angle  Run  Series  Results 

Table  4- 1  summarizes  the  results  of  the  runs  conducted  in  the  missile 
approach  angle  run  series.  The  initial  missile  approach  angle  is  found  in  column 
one.  For  each  approach  angle,  a  set  of  volume  plots  is  generated  (Appendix  A).  The 
best  cross-section  prediction  ( best  is  based  on  both  the  final  cross-section 
predicted,  and  the  true  missile  trajectory  being  within  the  filter  estimated  cross- 
section)  for  each  approach  angle  is  presented  in  column  six.  Each  volume  cross- 
section  prediction  is  made  by  removing  all  filter  measurements  and  allowing  the 
filter  covariance  estimates  to  diverge  (hereafter  referred  to  as  filter  divergence) 
for  one  second.  Columns  two  through  six  each  possess  two  sets  of  numbers.  The 
first  number  represents  the  filter  statistics  at  the  time  of  filter  covariance  release, 
and  the  second  set,  after  the  filter  is  allowed  to  diverge  for  one  second.  (For 
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example  row  two  of  Tabic  4-1  reads:  The  filter  took  measurements  until  t  =  5.0 
seconds  and  then  made  a  volume  prediction  for  t  =  (>.()  seconds.  At  5.0  seconds  the 
filter  estimated  range  is  900  (m)  with  an  estimated  rms  error  of  0.148  (m)  and  the 
estimated  LOS  rms  error  is  10.8  (urad).  After  1  second  of  propagation  into  the 
future,  the  range  is  estimated  as  200  (m)  with  a  filter  estimated  16.6  (m)  rms  error, 
and  with  a  filter  estimated  LOS  rms  error  of  4.18  (mrad).)  The  volume  cross-section 

Table  4-1  Summary  of  Missile  Approach  Angle  Run  Series  Results 


Approach 

Angle 

(radians) 

Time 

(seconds) 

A  Range 
(meters) 

A  Angle 
(mrad) 

Range 

(meters) 

HUH 

0.0 

BBS 

.152/12.6 

.0232  /  2.36 

1500/1050 

.0053/31.2 

jt/4 

5.0 -6.0 

.148/16.6 

.0168/4.18 

900  /  200 

.0022/13.9 

it/2 

4.0 -5.0 

.150/8.67 

.01 59/. 530 

1800/1000 

.0043  /  4.59 

3ut/4 

missile  failed  to  track 

X 

warn 

.150/11.7 

.01 57/. 982 

1700/700 

.0070  /  8.76 

Sk/4 

7.0-8.0 

.148/2.95 

.0171  /  2.72 

1900  /  800 

.0048/31.5 

3k/2 

BUHBI 

.149  / 19.0 

.0158  /  .368 

2300/1150 

.126/8.04 

7n/ 4 

missile  failed  to  track 

estimate  provided  in  column  six  is  calculated  based  on  Equation  (1-1)  with  two 
exceptions  -  a  small  angle  approximation  is  used  (since  the  largest  angle  in  Table  4- 
1  is  4  mrad  and  the  small  angle  approximation  greatly  simplifies  calculations),  and 
the  cross-section  values  are  based  on  lo  values  instead  of  the  3o  values  (again  just 
to  simplify  calculations)  used  in  Equation  (1-1).  Equation  4-1  presents  the  small 
angle  approximation  form  used  in  calculating  the  cross-sections: 

Volume  Cross-Section  =  (dSR)  (5K30)  (4-1) 
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where 


'Jt  is  the  filter  estimated  range  (col  5), 

diH  is  the  filter  estimated  range  rms  error  (col  £), 

d®  is  the  filter  estimated  angle  error  (col  4). 

Two  missile  approach  angle  run  sets,  the  n/4  and  5*/4  geometries,  are 
selected  for  additional  study  to  determine  possible  relationships  between  predicted 
volume  cross-section  and  missile  range.  The  two  missile  approach  angles  are 
selected  from  the  eight  missile  approach  angles  available  because  these 
simulations  have  longer  range  and  longer  simulation  run  times  than  the  other  six 
missile  approach  angles,  thus  providing  more  possible  data  points  for  study.  For 
each  of  these  two  approach  angles,  the  filter  is  run  for  a  specified  amount  of  time 
and  is  then  allowed  to  diverge  for  one  second  generating  a  cross-section 
prediction.  The  filter  run-time  (before  divergence)  starts  at  1.0  seconds  for  the 
first  run,  and  is  increased  by  0.5  seconds  for  each  additional  run.  The  filter  is 
propagated  to  each  specified  time  using  the  same  initial  conditions  and  random 
number  seeds,  and  is  then  diverged  for  one  second  to  generate  a  volume  cross- 
section  prediction,  it  is  thought  that  a  relationship  between  estimated  volume 
cross-section  and  range  can  be  established.  Because  the  volume  cross-section 
given  in  Equation  (4-1)  is  linear  versus  range,  a  "tilted-bathtub"  pattern  similar  to 
that  shown  in  Figure  4-1  is  expected.  The  "tilted-bathtub"  pattern  has  three 
distinct  parts,  the  initial  slope  which  represents  the  period  before  the  filter 
actually  converges  to  the  optimal  missile  states,  the  linear  region  where  the  filter 
covariance  estimates  are  assumed  to  be  stable  and  slowly  changing,  and  the  cross- 
section  estimate  decreases  inversely  with  range,  and  the  final  tail  where  the  range 
approaches  zero  causing  the  filter  to  go  unstable.  When  the  missile-to-aircraft 
range  approaches  zero,  small  changes  in  the  range  cause  very  large  changes  in 
the  LOS  and  LOS-rate.  This  causes  the  simulation  dynamics  to  go  numerically 
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Range 


0.0 


Figure  4-1  Expected  "Bathtub"  Pattern  for  Cross-Section  vs  Range 

unstable.  The  tracking  missile  seeker  experiences  the  same  problem.  To 
compensate  for  this  effect,  missiles  are  preprogrammed  with  a  special  flight 
trajectory  which  is  activated  when  the  missile-to-aircraft  tracking  errors  exceed  a 
prescribed  limit.  The  activation  of  the  flight  trajectory  causes  the  missile  to  enter 
the  "endgame  scenario"  in  which  it  no  longer  uses  proportional  navigation  as  a 
guidance  law.  In  the  filter,  this  stage  is  marked  by  a  sudden  rapid  growth  in  the 
filter  covariances  followed  by  an  integration  failure  within  MSOFE.  These  effects, 
which  ultimately  determined  the  scenario  time  length,  can  be  seen  in  Figure  A-86 
Appendix  A. 

Figure  4-2  (a  &  b)  shows  a  series  of  volume  cross-section  predictions  for  the 
jc/4  geometry.  Figure  4-2a  shows  the  volume  cross-section  at  filter  release  as  a 
function  of  range.  Figure  4-2  b  shows  the  volume  cross-section  after  one  second  of 
filter  divergence.  Note  that  the  subdivisions  on  the  range  axis  do  not  correlate 
from  Figure  4-2a  to  Figure  4-2  b.  The  range  axis  is  not  scaled  relative  to  range,  but 
rather  relative  to  time.  The  range  values  shown  are  taken  at  0.5  second  intervals 
beginning  at  1.0  seconds  in  Figure  4-2a.  Since  Figure  4-2b  represents  a  1.0  second 


4-4 


prediction,  its  range  axis  is  shitted  two  range  subdivisions  (or  one  second  in  time) 
to  the  right  For  instance  the  initial  cross-section  estimate  at  range  =  42(H)  (Figure 
4-2a)  corresponds  to  the  diverged  cross-section  estimate  at  range  =  4000  (Figure  4- 
2b)).  Figure  4-4  (a  &  b)  is  similarly  aligned.  Figures  4-2a  and  4- 3a  show  that  the 
expected  relationship  between  cross-section  estimation  and  range  is  achieved  at 
the  initial  filter  release.  However,  as  seen  in  Figures  4-2b  and  4-3b,  that  same 
relationship  does  not  exist  after  divergence.  Figure  4-3a  also  contains  two  outlier 
points.  The  mid-range  point  is  caused  by  a  brief  filter  tracking  error  and  is  not 
significant.  The  end-range  point  is  caused  by  the  "zero  range"  divergence  which 
is  now  starting  to  occur,  fhe  expected  "zero  range"  divergence  does  infact  occur 
but  can  not  be  shown  in  the  figures  because  when  it  occurs,  the  filter  divergence 
is  so  severe  it  crashes  the  computer  simulations  (usually  in  less  than  0.1  seconds). 
(For  an  example  of  range  induced  filter  divergence  see  Appendix  A  Figure  A-86.) 
In  the  first  section  of  the  "  bathtub"  curve,  the  filter  does  not  estimate  larger 
volumes,  but  instead  tends  to  lose  track  of  the  true  missile  trajectory  during 
divergence,  resulting  in  bad  cross-section  predictions. 


Figure  4-2a  Initial  Cross-Section  vs.  Range  (pi/4) 
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Figure  4-2b  Cross-Section  vs  Range  After  One  Second 
Divergence  (pi/4) 


Figure  4- 3  a  Initial  Cross-Section  vs  Range  (Spi/4) 


4.2  The  Performance  Baselines 

The  n/4  missile  approach  angle  is  used  in  all  the  remaining  simulations 
conducted.  The  n/4  geometry  (reference  Figure  3-2  for  a  pictorial)  is  selected 
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Figure  4-3b  Cross-Section  vs  Range  Alter  one  Second 
Divergence  (SFi/4) 


Range 


because  the  dynamics  in  that  geometry  are  least  sensitive  to  system  effects  likely  to 
cause  simulation  failure  (such  as  the  previously  discussed  zero  range  effect).  All 
simulations  are  conducted  using  identical  initial  conditions  and  random  number 
seeds  with  the  exception  of  the  parameter  currently  under  study.  The  baseline 
conditions  are  established  as  indicated  in  Table  3-2. 

Filter  plots  showing  range,  range-rate  angle  and  angle-rate  errors  are 
presented  in  Figure  4-4  (a-d).  Figure  4-4  (e-f)  presents  a  three  second  filter 
divergence  for  the  low  noise  baseline  initiated  three  seconds  after  simulation  start. 
Ail  remaining  simulations  are  conducted  using  identical  filter  divergence  criteria. 

A  high  noise  performance  baseline  is  similarly  established  using  the  n/ 4 
geometry  and  conditions  listed  in  Table  3-5.  The  high  noise  baseline  filter  plots  are 
presented  in  Figure  4-5  (a-d),  with  the  three-second  filter  divergence  plots  in 
Figure  4-5  (e-f)  respectively.  The  intention  is  to  establish  a  comparative 
performance  range  between  the  low  noise  baseline  results  and  the  high  noise 
baseline  results.  The  primary  figures  of  merit  for  the  comparison  are  the  filter 
estimated  rms  errors  after  one  second  and  three  seconds  of  filter  divergence. 
Throughout  this  analysis  true  errors  are  only  considered  to  the  extent  that  they 
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Figure  4-4  (a  and  b) 
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Figure  4-4  (c  and  d) 
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Figure  4-4  (e  and  f) 
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Figure  4-5  (a  and  b) 
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are  sufficiently  ( by  at  least  one  standard  deviation)  outside  the  filter's  estimated 
rms  errors  indicating  an  invalid  cross-section  prediction.  The  use  of  filter 
estimated  statistics  is  as  an  analysis  technique  is  used  throughout  this  thesis  as  a 
method  of  analyzing  the  filter  predicted  volume.  It  is  not  the  intent  of  this  thesis  to 
study  issues  particular  to  optimal  extended  Kalman  filter  design.  More  will  be  said 
about  the  results  of  this  assumption  in  Chapter  V.  The  high  noise  baseline 
divergence  plots  only  go  out  to  five  seconds.  The  three  second  divergence  goes 
unstable  shortly  after  two  seconds  of  propagation,  so  the  three  second  worst  case 
performance  can  not  be  established.  The  filter  estimated  rms  error  for  the 
baseline  simulations  are  summarized  in  Table  4-2. 

fable  4-2  Baseline  Performance  Results 


Baseline 

Noise 

One  Second 
rms  Range 
Error 

Divergence 
rms  Angle 
Error 

Three  Secorn 
rms  Range 
Error 

i  Divergence 
rms  Angle 
Error 

low 

6.336  (m) 

537.5  (urad) 

67.61 

5.33  mrad 

high 

10.4  (m) 

1.11  (mrad) 

failed  n/a 

failed  n/a 

4.3  The  Ladar  Run  Series  Results 

The  parameters  studied  in  the  ladar  series  run  sets  are  summarized  in  Tables 
3-2  and  3-3.  The  results  are  summarized  in  Tables  4-3  (a-e).  All  values  presented  in 
the  first  column  are  the  rms  noise  values  used  in  simulation.  The  values  marked 
with  an  "L"  indicate  that  the  truth  model  (only)  uses  the  values  given  in  the  table 
while  the  filter  remains  at  baseline  conditions.  The  values  marked  with  an  "H" 
indicate  that  the  filter  model  (only)  uses  the  values  given  in  the  table  while  the 
truth  model  remains  at  baseline.  The  unmarked  column  one  entries  indicate  that 
both  the  truth  model  and  the  filter  model  use  these  values.  Columns  two  through 
seven  give  the  filter  estimated  rms  errors  at  initial  filter  release,  after  one  second 
of  filter  divergence  and  after  three  seconds  of  filter  divergence. 
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The  result  of  the  individual  ladar  measurement  study  indicates  that  no  one 
ladar  measurement  dominates  the  overall  volume  prediction  at  either  one  second  or 
three  seconds  of  filter  divergence.  Similarly  dropping  the  ladar  pulse-rate  to  SO  Hz 
also  does  not  significantly  affect  the  figures  of  merit,  fable  4-3  shows  that,  with 
the  exception  of  two  failed  angle  simulations,  the  final  predicted  volume  after 
three  seconds  of  divergence  never  consistently  exceeds  50%  of  the  value 
established  by  the  baseline,  and  in  some  cases  even  bettered  the  established 
performance  baseline.  The  results  suggest  that  the  variation  in  predicted  volume 
as  a  function  of  ladar  parameter  variation  is  not  large  enough  to  be  statistically 
significant  using  only  filter  estimated  statistics.  Monte-Carlo  simulations  (which 

Table  4- 5a  Ladar  Range  Results 


RMS 

Noise 

Value 

Initial  Filter  Release 

One  Second 

Three  Seconds 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

1.0  m 

0.129 

163 

6.10 

461 

59.4 

3690 

2.0  m 

0.193 

17.7 

735 

802 

91.85 

8180 

■9 

0327 

163 

6.66 

622 

72.16 

5850 

1.0  mL 

.129 

173 

6.91 

713 

77.9 

6380 

2.0  mL 

.193 

16.9 

6.11 

546 

613 

4520 

4.0  mL 

327 

17.6 

7.18 

753 

801 

7390 

1.0  mH 

.0906 

17.5 

720 

756 

853 

6860 

2.0  mH 

.091 

17.2 

6.12 

584 

62.1 

4690 

4.0  mH 

.091 

17.0 

633 

528 

164 

4320 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians. 
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Table  4-3b  ladar  Range-Rate  Results 


RMS 

Noise 

Value 

Initial  Filter  Release 

One  Second 

Three  Seconds 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Lrror 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

2.0  m/s 

0.133 

16.9 

7.63 

701 

87.2 

7540 

4.0  m/s 

0.165 

16.4 

7.49 

476 

65.16 

3930 

8.0  m/s 

.199 

16.8 

7.83 

447 

64.26 

3360 

2.0  m/s  L 

.125 

17.1 

7.13 

620 

72.2 

mm 

4.0  m/s  L 

.166 

17.5 

mm 

855 

110 

10400 

8.0  m/s  L 

.199 

17.1 

8.15 

531 

68.7 

4390 

2.0  m/s  H 

.0906 

17.4 

634 

701 

80.1 

6460 

4.0  m/s  H 

.091 

17.1 

7.16 

750 

98.7 

7670 

8.0  m/s  H 

.0905 

13.6 

933 

37 7 

130 

2440 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians. 


Table  4-3 c  Ladar  Angle  Results 


Initial  Filter  Release 

One  Second 

Three  Seconds 

RMS 

RMS 

RMS 

RMS 

RMS 

RMS 

RMS 

Noise 

Range 

Angle 

Range 

Angle 

Range 

Angle 

Value 

Error 

Error 

Error 

Error 

Error 

Error 

80  mad 

.0906 

28.2 

639 

630 

6731 

5400 

KSSH99 

.0906 

453 

539 

458 

49.5 

2990 

250  urad 

divergence  failed 

80  mad  L 

.0906 

283 

634 

746 

842 

7460 

160  mad  L 

.0905 

47.7 

7.10 

845 

923 

7920 

divergence  failed 

80  urad  H 

.0906 

17.4 

638 

7 07 

81.1 

7010 

160  urad  H 

.091 

173 

632 

691 

832 

6430 

250  urad  H 

0902 

163 

522 

374 

48.6 

6430 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians. 


Table  4-3d  Ladar  Angle-Rate  Results 


RMS 

Noise 

Value 

Initial  Filter  Release 

One  Second 

Three  Seconds 

BsSSSH 

RMS 

Angle 

Error 

RMS 

Range 

Frror 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

4.0  mr/s 

.0907 

17.8 

6.78 

681 

743 

6150 

8.0  mr/s 

.0907 

17.5 

621 

595 

643 

5020 

4.0  mr/s  L 

.091 

13.7 

6.40 

562 

65.0 

8.0  mr/s  L 

.091 

17.6 

635 

warn 

66.4 

5170 

4.0  mr/s  H 

.0905 

16.8 

627 

611 

70.3 

5470 

8.0  mr/s  H 

.0905 

17.1 

wmm 

715 

863 

7380 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians. 


Table  4-3e  Ladar  Pulse-Rate  Results 


RMS 

Noise 

Value 

Initial  Filter  Release 

One  Second 

:m  m  u  ,i — j,  ,i  ,  '  ■„  ,m  i 

Three  Seconds 

EH 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

60  Hz-2 

.180 

30.1 

733 

761 

77.9 

6240 

60  Hz-3 

350 

53.1 

820 

864 

873 

6990 

60  Hz-4 

.720 

79.1 

9.81 

966 

88.1 

7320 

30  Hz-1 

.125 

213 

736 

750 

813 

7250 

30  Hz-2 

.247 

36.0 

8.11 

616 

723 

5070 

30  Hz-3 

.494 

70.1 

9.75 

1008 

423 

8380 

30  Hz-4 

.986 

99.5 

10.6 

864 

793 

5620 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians. 


unfortunately  are  not  used  in  this  thesis  study)  are  required  to  establish  any 
meaningful  conclusions  on  the  relationships  between  predicted  volume  size  and 
rms  ladar  noise. 
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There  are  two  notable  exceptions  to  these  results.  The  simulation  dynamics 
are  more  sensitive  to  rms  angle  noise,  than  to  rms  range,  range-rate  and  angle- 
rate  noise  as  is  seen  from  Table  -Me.  The  angle  noise  sensitivity  does  not  appear  in 
the  volume  predictions,  but  rather  within  the  system  dynamics  themselves.  At 
higher  rms  angle  noise  the  filter  becomes  unstable. 

The  second  notable  exception  is  found  by  close  analysis  of  the  results 
presented  in  Table  4-4e.  For  each  ladar  pulse-rate,  four  cases  are  tested  (60  Hz-1  is 
the  low  noise  baseline  presented  in  Table  4-2)  using  the  rms  measurement  noise 
sets  from  Table  (3-2).  At  the  initial  filter  release  time,  there  is  a  strong 
relationship  between  estimated  volume  and  noise  strength.  A  similar  relationship 
exists  between  the  filter  estimated  volume  and  ladar  pulse-rate  at  filter  release.  For 
the  60  Hz  pulse-rate,  the  relationship  between  volume  and  ladar  noise,  although 
still  present,  is  less  prevalent  for  one  second  of  filter  divergence  and  even  more  so 
after  three  seconds  of  filter  divergence.  For  the  30  Hz  pulse-rate,  there  is  a 
marginal  relationship  between  ladar  noise  and  filter  estimated  range  error  after 
one  second  of  filter  divergence.  There  is  no  evidence  of  relationship  between 
filter  estimated  range  error  and  ladar  noise  for  30  Hz  after  three  seconds  of 
divergence.  There  is  no  evidence  of  relationship  between  filter  estimated  angle 
error  and  ladar  noise  for  either  divergence  at  30  Hz.  Further  discussion  of  Table  4- 
3e  results  is  deferred  to  Chapter  V. 

4.4  The  Aircraft  Run  Series  Results 

Simulations  in  the  aircraft  run  series  are  conducted  using  the  criteria 
established  in  Section  3.3  and  Table  3-4.  For  consistency  a  re/4  missile  approach 
angle  is  used,  as  discussed  in  Section  4-2.  Throughout  the  aircraft  run  series,  no 
aircraft  acceleration  measurements  are  assumed  to  be  available  from  the 
navigation  system.  All  other  parameters  (except  the  one  under  study  in  the 
particular  run  set)  are  held  at  the  baseline  values.  The  aircraft  run  series 
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investigates  the  effects  ot  assumed  INS  measurement  noise,  assumed  pilot  error 
noise,  and  aircraft  turn  rate  on  the  predicted  volume  cross-section. 

The  low  noise  baseline  results  presented  in  Table  4-2  assume  a  (0.  lg)2/sec 
pilot  induced  aircraft  trajectory’  error.  The  baseline  value  is  taken  from  the 
development  in  Chapter  11.  Since  the  baseline  value  is  based  on  an  ad  hoc  estimated 
noise  strength  (see  the  discussion  on  pilot  induced  error  in  Section  2.4),  this 
parametric  analysis  is  performed  to  assess  the  impact  of  assumed  pilot  induced 
error  on  the  predicted  cross-section.  Table  4-4  presents  the  results  obtained  as  the 
assumed  pilot  induced  error  is  increased.  From  the  Table  4-4  results,  two 
observations  are  made.  The  filter's  cross-section  estimate  at  filter  release  appears 
to  be  constant,  indicating  that  the  velocity  measurements  provided  by  the 
navigation  subsystem  are  capable  of  filtering  out  the  pilot  induced  noise  strength 
as  long  as  measurements  are  available.  For  the  one-second  and  three-second  filter 
divergences,  there  is  no  obvious  relationship  between  pilot  induced  error  and 
estimated  cross-section,  indicating  the  need  for  Monte-Carlo  analysis. 

Table  4-5  presents  the  results  of  the  INS  run  set  and  Table  4-6  presents  the 
results  of  the  aircraft  turn  rate  run  set.  The  results  are  consistent  with  those 


Table  4-4  Pilot  Induced  Error  Results 


Mean 

Squared 

Noise 

Strength 

Initial  Filter  Release 

One.  Second 

Three  Seconds 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

(0.2g)2/s 

.0906 

17.1 

5.88 

934 

583 

3940 

(0.3g)2/s 

.0906 

17.2 

7.17 

761 

963 

7880 

(0.4g)2/s 

.0906 

16.8 

626 

606 

71.6 

5490 

{OSS)2/* 

.0906 

17.4 

627 

576 

67.1 

4610 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians. 
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observed  in  the  pilot  induced  error  run  set  (i.e.,  no  relationship  between 
parametric  variations  and  predicted  volume  cross-section  at  either  one  or  three 
seconds  of  divergence).  For  purposes  of  comparison,  the  baseline  values  presented 
in  Table  4-2  include  an  assumed  0.1  m/s  velocity  error  (with  acceleration 
measurements)  and  a  4.0  g  aircraft  turn.  Table  4-5  results  for  the  cross-section  at 
initial  filter  release  shows  that  filter  estimated  range  and  angle  errors  do  not 
change  with  an  increase  in  assumed  aircraft  velocity  measurements.  The  same 
results  from  Table  4-5  indicate  that  it  may  be  possible  to  eliminate  acceleration 
measurements  from  the  filter  model  altogether,  as  performance  without 
acceleration  measurements  is  only  marginally  above  baseline  performance  (at 
initial  filter  release).  Table  4-6  results  at  initial  filter  release,  show  that  aircraft 
turn  rate  does  affect  initial  cross-section  estimates.  Neither  table  shows  a 
relationship  between  the  parameter  under  study  and  the  filter  estimated  cross- 
section  at  either  one  or  three  seconds  of  divergence,  again  indicating  the  need  for 
Monte-Carlo  simulation. 


4.5  The  PRONAV  Parameter  Pseudo-Noise 

The  results  of  the  parametric  studies  conducted  to  this  point  yield  no 
conclusive  results  concerning  which  model  parameter  has  the  largest  impact  on 


Table  4-5  INS  Velocity  Measurement  Error  Results 


RMS 

Velocity 

Error 

Initial  Filter  Release 

One  Second 

Three  Seconds 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

RMS 

Range 

Error 

RMS 

Angle 

Error 

0.5  m/s 

.0906 

175 

7.40 

798 

992 

8480 

1.0  m/s 

.0906 

175 

656 

689 

765 

6320 

5.0  m/s 

.0906 

17.4 

6.76 

677 

762 

6190 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians 
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Table  4-6  Aircraft  turn  Kate  Results 


Aircraft 
Turn  Rate 

Initial  Filter  Release 

One  Second 

Three  Seconds 

RMS 

Range 

Frror 

RMS 

Angle 

Frror 

RMS 

Range 

Frror 

RMS 

Angle 

Frror 

RMS 

Range 

Frror 

RMS 

Angle 

Frror 

1.0  g  turn 

.0906 

16.4 

6.72 

-m 

70.3 

4230 

2.0  g  turn 

.0906 

15.7 

4.82 

200 

38.6 

7880 

4.0  g  turn 

.0906 

17.1 

626 

606 

71.6 

5490 

5.0  g  turn 

.0910 

18.2 

6.27 

576 

67.1 

4610 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians 


the  divergence  of  the  volume  cross-section.  This  prompts  further  simulation  to 
research  exactly  what  model  parameter  does  have  a  dominant  impact  on  the  filter 
divergence.  The  ladar  measurement  noise,  the  INS  noise,  the  pilot  induced  error, 
the  aircraft  maneuver,  and  the  scenario  conditions  have  all  been  subjected  to 
parametric  analysis  and  none  have  shown  a  strong  relationship  between  the 
parametric  variations  and  the  filter  predicted  cross-section  after  either  one  or 
three  seconds  of  filter  divergence.  Only  three  possible  candidates  remain,  the 
proportional  navigation  (Pro-Nav)  constant  (n)  pseudo-noise,  the  filter's  surface- 
to-mass  ratio  (S/M)  pseudo-noise,  and  the  lateral  acceleration  tuning  noise.  The 
Pro-Nav  constant's  pseudo-noise  turns  out  to  be  the  driver.  (Like  the  parameters 
studied  earlier,  the  surface-to-mass  ratio  and  lateral  acceleration  noise  strengths 
have  no  discemabie  effect  on  the  predicted  volume  cross-section.)  The  Pro-Nav 
constant  is  given  an  initial  diffusion  strength  Q(t)  of  (.5)^/sec.  This  Q{t)  value 
allows  the  filter  to  occasionally  lose  track  of  the  system  dynamics,  particularly  in 
harsher  dynamic  scenarios  where  the  missile  require  larger  lateral  accelerations. 
The  Qlt)  value  is  increased  to  (2)Vsec  to  facilitate  better  lateral  acceleration 
tracking  (which  it  does).  However,  on  filter  release  the  Pro-Nav  constant  diffusion 
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strength  is  now  dominating  the  covariance  spread  (and  hence  the  predicted 
volume  cross-section),  lo  investigate,  the  aircrait  turn  rate  set  is  resimulated 
using  a  (i(t)  of  (,5)2/sec  for  the  Pro-Nav  constant.  Ihe  results  of  the  new  run  set 
are  presented  in  Table  4-7.  A  second  run  set  is  simulated  using  the  baseline 
conditions  and  varying  the  Pro-Nav  constant  noise  diffusion  strength  for  direct 
comparison.  Ihese  results  are  presented  in  Table  4-8. 

Table  4-7  Aircraft  Turn  Kate  Results  Using  New  Pro-Nav  Constant  Qj(t) 


Aircraft 
Turn  Rate 

Initial  filter  Release 

One  Second 

Three  Seconds 

RMS 

Range 

lirror 

RMS 

Angle 

lirror 

RMS 

Range 

lirror 

RMS 

Angle 

lirror 

RMS 

Range 

lirror 

RMS 

Angle 

frror 

0.1  g  turn 

.0906 

13.7 

6.70 

439 

69.1 

1720 

1.0  g  turn 

.0906 

13.4 

528 

182 

50.4 

1100 

2.0  g  turn 

.0906 

13.5 

4.82 

200 

38.6 

1320 

3.0  g  turn 

.0906 

13.8 

171 

27.5 

1090 

4.0  g  turn 

divergence  failed  0.5  seconds  after  release 

5.0  g  turn 

divergence  failed  0.26  seconds  after  release 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians 


Table  4-8  n  Parameter  Q(t)  vs  Baseline  Comparison 


Oft) 

Initial  Filter  Release 

One  Second 

Three  Seconds 

Range 

Angle 

Range 

Angle 

Range 

Angle 

baseline 

.0906 

15.4 

634 

538 

67.6 

5330 

(1.0)2/sec 

.0905 

153 

521 

370 

45.4 

3020 

(0.5 )2/sec 

.0906 

133 

4.09 

171 

273 

1320 

All  range  values  are  given  in  units  of  meters. 

All  angle  values  are  given  in  units  of  micro-radians 
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The  results  in  Table  4-8  show  a  definite  correlation  between  the  Pro-Nav 
constant  diffusion  strength  and  the  filter  divergence.  Although  the  divergence 
correlation  is  strong  for  both  range  and  angle  divergence,  there  is  stronger 
divergence  in  the  angle,  fable  4-7  results  verify  that  the  filter's  estimated  cross- 
section  is  significantly  smaller  at  all  divergences  when  using  a  lower  Pro-Nav 
constant  Qlt).  Table  4-7  also  verifies  that  the  lower  Q{t)  value  does  not  allow 
adequate  tracking  when  the  scenario  dynamics  require  higher  missile  lateral 
acceleration  as  seen  from  the  failed  divergences  under  aircraft  4.0  g  and  5.0  g 
turns. 

4.6  Summary 

This  chapter  presents  the  simulation  results  obtained  from  the  kill  zone 
volume  parametric  analysis  described  in  Chapter  III.  In  general  the  initial  volume 
cross-section  at  filter  release  varied  as  expected  with  the  parameter  under  study. 
However,  the  volume  cross-section  after  three  seconds  of  divergence  is,  in  general, 
independent  of  the  parameter  under  study.  New  simulations  are  conducted  to 
investigate  why  the  filter  predicted  three  second  cross-section  is  independent  of 
ladar  measurement  noise,  assumed  INS  rms  errors,  pilot  error  and  aircraft  turn 
rate.  It  is  found  that  the  tuning  pseudo-noise  used  on  the  filter's  estimation  of  the 
proportional  navigation  parameter  (n)  is  dominating  the  filter  divergence.  The 
other  parameters  effects  can  not  be  adequately  ascertained  without  Monte-Carlo 
simulation. 

Chapter  V  summarizes  the  research  accomplished  in  this  thesis  study.  The 
work  presented  in  Chapters  I  to  IV  is  briefly  reviewed,  the  main  conclusions  of  the 
thesis  study  are  presented,  and  recommendations  for  continued  research  are  made 
based  on  these  conclusions.  In  particular,  certain  results  presented  in  this  chapter 
are  reviewed  in  light  of  the  Pro-Nav  constant  induced  filter  divergence,  and  in 
light  of  the  obvious  need  for  Monte-Carlo  simulation. 
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V  Conclusions  and  Recommendations 


5.0  Introduction 

The  conclusions  presented  in  this  chapter  are  based  on  the  results  presented 
in  Chapter  IV.  The  assumptions  and  analysis  techniques  used  to  scope  this  thesis 
study  are  summarized.  The  conclusions  focus  on  the  impact  of  the  Pro-Nav 
parameter's  noise  diffusion  and  its  effect  on  the  other  simulation  results.  Research 
recommendations  are  made  based  on  the  conclusions  presented  in  this  chapter,  the 
assumptions  outlined  in  Chapter  1  and  the  model  limitations  from  Chapter  II.  The 
chapter  ends  with  an  overall  thesis  summary. 

5.1  Review  of  Thesis  Assumptions 

Proper  interpretation  of  the  results  presented  in  Uiapter  IV  first  requires  a 
review  of  the  assumptions  made  in  Chapters  I  and  II  which  scope  the  thesis  and 
define  the  limitations  of  the  simulation  model.  Therefore,  as  a  first  step  to 
presenting  conclusions,  the  modelling  assumptions  on  which  those  conclusions  are 
based,  are  now  presented. 

A  benign  and  easily  tracked  aircraft  trajectory  is  chosen  to  simplify  model 
dynamics.  The  simple  trajectory  also  reduces  required  missile  lateral  acceleration, 
obstensibly  making  the  filter's  tracking  and  predicting  functions  easier.  From  a 
mission  standpoint  flying  a  benign  trajectory  may  not  be  desirable  or  even 
possible. 

The  analyses  performed  throughout  this  thesis  use  filter  computed  statistics 
as  the  benchmark  of  performance.  This  may  be  true  of  an  optimally  designed  and 
tuned  filter,  which  has  been  thoroughly  tested  with  Monte-Carlo  analysis,  but  no 
such  design  and  analysis  has  been  performed  as  part  of  this  thesis.  This  thesis 
looks  past  optimal  filter  design  issues  to  get  a  first  glimpse  at  how  the  filter 
computes  and  propagates  the  predicted  volume  cross-section.  This  thesis  study  is 
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intended  as  only  the  first  in  a  series  of  studies  on  the  complex  subject  of  kinetic  kill 
zone  prediction,  as  such  specific  design  issues  have  been  surpassed  to  get  a  big 
picture  look  at  the  whole  volume  prediction  problem. 

The  missile  model  used  in  this  thesis  assumes  a  post  burnout  coast-to- 
intercept  missile  flight  with  no  endgame  guidance.  As  such  the  missile  simulations 
conducted  are  limited  to  those  in  which  a  kineticaliy  limited  missile  is  capable  of 
adequately  tracking  the  target.  The  envelope  of  possible  engagement  scenarios 
can  be  significantly  expanded  by  using  a  thrusting  missile  model  in  future 
analyses.  The  missile  endgame  guidance  should  also  be  considered  particularly 
when  using  a  thrusting  missile  model.  Most  missiles  have  lateral  acceleration 
limits  between  24  g  and  40  g.  The  non-thrusting  missile  is  not  kineticaliy  capable 
of  exceeding  about  14  g  in  lateral  acceleration  before  bleeding  off  airspeed  and 
stalling.  A  thrusting  missile  may  be  capable  of  exceeding  this  limit  and  therefore 
needs  to  consider  the  addition  of  an  endgame  trajectory  when  specific  acceleration 
limits  are  exceeded. 

The  initial  ladar  rms  measurement  errors  developed  in  Chapter  n  are  based 
on  a  benign  atmospheric  turbulence  model.  A  good  jet-wash  turbulence  model  has 
not  yet  been  developed.  More  adequate  measurement  statistics  based  on  the  active 
turbulence  caused  by  jet-wash  and  missile  flight  needs  to  be  added  to  the  ladar 
model.  This  thesis  overcomes  the  model  limitation  by  treating  the  assumed  ladar 
measurement  errors  as  parameters  of  study. 

The  simulations  are  limited  to  a  constant  altitude  plane  throughout  this 
study.  The  baro-altimeter  (vertical  position)  measurements  used  in  this  thesis,  are 
useful  within  the  context  of  a  constant  altitude  simulation,  but  are  not  realistic  if 
the  missile  is  not  confined  to  the  aircraft's  altitude  plane.  A  three  dimensional 
simulation  definitely  addresses  the  issue  of  observing  missile  altitude.  Should 
however,  a  two  dimensional  model  be  selected  for  continued  analysis,  the  altimeter 
measurements  should  be  removed  from  the  model. 
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During  simulation  the  missile  truth  model  used  a  Pro-Nav  constant  of  4.5. 

The  filter  estimate  of  the  Pro-Nav  constant  consistently  settled  to  -4.5.  A  thorough 
search  was  made  to  hunt  for  a  possible  sign  error  in  the  code  to  explain  the 
negative  filter  estimate,  but  came  up  empty.  The  simulation  results  indicate  that 
the  missile  does  track  the  aircraft  in  both  the  filter  and  truth  models  (see 
trajectory  plots  in  Appendix  A).  All  tests  which  attempted  to  switch  the  sign  bias  of 
the  filter  model  led  to  almost  immediate  simulation  because  the  filter  initially  tried 
to  laterally  accelerate  in  the  wrong  direction.  Thus  the  negative  bias  for  the  filter 
estimate  of  the  Pro-Nav  gain  remains  unexplained.  It  is  mentioned  here  as  a 
warning  to  future  researchers  who  may  be  interested  in  using  the  FORTRAN  code 
developed  during  this  study. 

5.2  Conclusions 

Discussion  now  returns  to  the  ladar  run  series  results  presented  in  Table  4- 
3c  and  Table  4-3e.  It  is  believed  that  the  angle  measurement’s  250  i*rad  rms  error 
runs  fails  during  filter  divergence  because  the  Pro-Nav  constant's  noise  diffusion 
has  its  largest  impact  on  angle  divergence.  Spreading  the  ladar  angle 
measurements  while  maintaining  tight  ladar  measurements  of  range,  range-rate 
and  angle-rate,  causes  an  imbalance  in  the  filter  covariance  estimates  which  in- 
tum  causes  filter  instability  during  propagation  without  measurement. 

The  data  from  Table  4-3e  is  reconsidered  under  scrutiny  of  the  Pro-Nav 
constant's  noise  diffusion.  The  data  is  presented  in  a  new  format  in  Table  5-1. 
Column  one  of  the  table  lists  the  runs  being  compared.  The  notation  60  Hz  1-2 
indicates  that  the  results  from  the  60  Hz  case  1  run  are  being  compared  to  the 
results  from  the  60  Hz  case  2  run  (Table  4-3e).  The  notation  30-60  (1)  indicates  that 
the  30  Hz  case  1  run  is  being  compared  to  the  60  Hz  case  1  run.  The  results  indicate 
that  at  filter  release  time  the  filter  estimated  rms  range  and  angle  errors  basically 
follow  the  ladar's  assumed  rms  measurement  errors.  In  each  case  as  the  assumed 
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rms  error  is  doubled,  the  filter  estimated  range  error  doubles  and  the  filter 
estimated  angle  error  almost  doubles.  The  results  also  indicate  that  there  is  an 
approximate  30%  increase  in  filter  estimated  errors 


Table  5-1  Range  and  Angle  Divergence  Correlations 


Runs 

Initial  Filter  Release 

One  Second 

Three  Seconds 

Compared 

Range  A% 

Angle  A% 

Range  A% 

Angle  A% 

Range  A% 

Angle  A% 

60  Hz  1-2 

98.7% 

95.5% 

11.6% 

42.0% 

11.5% 

17.1% 

60  Hz  2-3 

94.4% 

76.4% 

11.9% 

13.5% 

11.2% 

12.0% 

60  Hz  3-4 

105.7% 

49.0% 

19.6% 

11.8% 

0.92% 

4.7% 

30  Hz  1-2 

97.6% 

69% 

3.2% 

no  corr 

no  corr 

no  corr 

30  Hz  2-3 

100% 

94.7% 

mm 

no  corr 

no  corr 

no  corr 

30  Hz  3-4 

99.6% 

41.2% 

8.7% 

no  corr 

no  corr 

no  corr 

30-60(1 

9 

38.0% 

38.3% 

no  corr 

no  corr 

no  corr 

no  corr 

30-60(2 

1 

37.2% 

19.6% 

no  corr 

no  corr 

no  corr 

no  corr 

30-60(3 

1 

41.1% 

32.0% 

no  corr 

no  corr 

no  corr 

no  corr 

30-60  (4 

1 

36.9% 

25.8% 

no  corr 

no  corr 

no  corr 

no  corr 

when  the  pulse-rate  is  decreased  from  60  Hz  to  30  Hz.  After  one  second  of 
divergence,  however  the  percentage  change  in  filter  estimated  errors  have 
dropped  to  about  15%  for  60  Hz  and  10%  (range  only)  for  30  Hz  as  the  assumed  ladar 
measurement  errors  double.  The  Table  results  also  show  that  doubling  the  assumed 
ladar  measurement  accuracy  yields  less  improvement  in  volume  prediction  at  30  Hz 
than  at  60  Hz. 

The  conclusion  that  the  Pro-Nav  constant's  pseudo-noise  strength  drives  the 
filter  divergence  without  measurements,  supports  the  observations  made  in  Section 
4.1.  Figure  4-2b  and  Figure  4-3b  both  indicate  that,  at  filter  release,  a  relationship 
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between  estimated  range  and  estimated  volume  size  exists  as  expected.  The  same 
relationship  is  not  observed  alter  one  second  of  filter  divergence.  As  a  final  note, 
the  claim  is  not  being  made  that  no  relationship  exists  between  filter  estimated 
ladar  measurement  noise  and  filter  estimated  cross-section,  or  estimated  range  and 
estimated  cross-section,  but  only  that  the  relationships  are  small  compared  to  that 
observed  from  variations  in  the  Pro-Nav  constant's  diffusion  strength. 

Determining  the  impacts  of  either  the  assumed  rms  ladar  noise  or  the  filter 
estimated  range,  on  the  post  divergence  filter  estimated  cross-section,  requires 
Monte-Carlo  simulation  which  is  not  done  in  this  study. 

This  thesis  investigates  the  use  of  an  extended  Kalman  filter  (EKF)  to  predict 
the  kinetic  kill  zone.  The  EKF  approach  is  not  the  only  valid  approach  to  such  a 
complex  problem.  Although  the  initial  EKF  approach’s  results  obtained  in  this 
study  are  not  considered  favorable,  the  EKF  approach  cannot  be  ruled  out  as  invalid 
based  on  Section  5.1  assumptions  (i.e  no  Monte-Carlo  simulations  are  conducted,  and 
the  filter  used  is  not  optimally  designed  and  tuned  for  the  kinetic  kill  zone 
problem). 

5.3  Recommendations 

The  EKF  approach  to  the  kinetic  kill  zone  estimation  problem  represents 
only  one  of  many  possible  approaches.  Given  the  results  presented  and  their 
ensuing  conclusions,  further  research  is  necessary  to  determine  if  the  EKF 
approach  is  an  optimum  approach.  There  are  considerable  processing  and  time 
line  risks  associated  with  an  EKF  approach  only  a  few  of  which  are  presented  in 
Chapter  I.  Despite  the  risk,  the  EKF  approach  is  used  because  it  is  thought  that  an 
EKF  approach  offers  performance  advantages  over  other  possible  approaches. 
Given  the  resulting  dominance  of  the  filter's  Pro-Nav  constant  noise  strength,  it  is 
not  clear  whether  an  optimally  designed  EKF  will  offer  these  performance 
advantages.  In  light  of  this,  it  is  recommended  that  alternative  approaches  to  the 
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kinetic  kill  zone  estimation  problem  be  studied. 

There  are  three  things  which  are  relatively  certain  about  the  kinetic  kill 
zone  estimation  problem.  They  are  ( 1 )  at  any  given  time,  a  single  ladar 
measurement  represents  a  reasonably  accurate  measurement  of  current  missile 
location  and  velocity,  (2)  the  missile  is  guiding  itself  towards  the  target  using  some 
form  of  proportional  navigation,  and  (3)  the  missile  must  traverse  the  space 
between  the  missile's  current  position  and  the  aircraft  estimated  trajectory  in 
order  to  intercept  the  aircraft. 

A  possible  new  approach  could  perform  Monte-Carlo  analysis  on  the  missile 
flight  path  itself  in  order  to  establish  a  mean  and  estimated  rms  missile  flight 
trajectory.  The  aircraft  and  missile  starting  locations  remain  fixed,  while  missile 
dynamic  models  are  subjected  to  random  variations,  to  generate  statistics  on  the 
possible  missile  trajectory  to  intercept.  The  results  can  be  used  as  a  database  to  aid 
in  the  kill  zone  estimation  problem.  The  only  online  estimation  that  occurs,  is 
linear  projection  of  the  missile  location  based  on  the  ladar  measurements.  The 
processing  timeline  risk  goes  away  since  such  calculations  can  probably  be 
accomplished  within  a  few  milliseconds.  Maximum  use  is  then  made  of  the  ladar 
measurements  accuracy  since  system  reaction  can  be  almost  instantaneous.  This 
allows  maximum  latitude  in  selecting  an  optimum  intercept  time  and  location 
which  is  not  driven  the  by  stringent  processing  timelines  that  an  EKF  approach 
requires. 

Should  the  EKF  approach  to  the  kinetic  kill  zone  estimation  problem  be 
explored  further,  the  following  recommendations  are  made. 

(1)  A  need  exists  for  Monte-Carlo  simulation  to  accurately  assess  the  impacts 
of  parametric  variations  on  model  performance.  Specific  conclusions  about  the 
effects  of  rms  ladar  measurement  errors,  pilot  induced  trajectory  error,  aircraft 
turn-rate,  and  estimated  rms  INS  errors  on  the  final  estimated  cross-section  could 
not  be  made  based  solely  on  the  results  of  filter  computed  statistics. 
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(2)  The  possibility  of  (1-switching  the  Pro-Nav  constant’s  pseudo-noise 
strength  needs  study  to  determine  if  a  high  initial  diffusion  strength  can  be  used  to 
enhance  tracking  performance,  and  a  lower  diffusion  strength  used  for  volume 
estimation.  This  may  require  the  aircraft  to  perform  specific  or  restricted 
maneuvers  to  enhance  the  performance  of  the  CL-switched  filter  (which  inturn 
may  not  be  desirable  from  a  mission  standpoint  as  discussed  earlier). 

(3)  This  study  uses  a  post  burnout  coast- to- intercept  missile  scenario.  The 
following  missile  model  changes  are  recommended:  The  missile  should  be  initially 
thrusting,  with  the  possibility  of  going  into  coast;  The  simulated  missile  should 
have  an  endgame  guidance  capability  to  allow  more  complete  simulations.  The 
implementation  of  three  dimensional  simulations  can  be  accomplished  at  a  more 
advanced  stage. 

(4)  This  study  uses  a  benign  turbulence  model  and  accounts  for  the 
inadequacies  in  the  model  by  performing  a  parametric  study  on  assumed  ladar 
noises.  A  more  accurate  representation  of  aircraft  turbulence  should  be  included 
in  the  next  study.  In  the  event  that  an  accurate  jet  wash  turbulence  model  does  not 
exist,  measurements  need  to  be  made  to  determine  turbulence  induced 
measurement  errors  more  completely. 

(5)  Processing  time  requirements  may  be  a  potential  problem  when  using 
the  EKF  approach  to  kinetic  kill  zone  estimation,  especially  as  advanced  models  are 
developed.  The  timing  diagram  presented  in  Figure  1-3  requires  online  filter 
estimation  in  faster  than  real  time.  Given  the  complexity  of  the  model  used  in 
simulation  and  additional  complexities  which  would  become  part  of  an  actual 
system,  it  is  not  certain  that  faster  than  real  time  simulation  can  be  achieved. 
Alternatives  to  filter  processing  may  need  investigation  just  on  this  point  alone. 

(6)  Accurate  assessment  of  time-to-go  requires  an  ability  to  estimate  missile 
endgame  guidance.  Any  time-to-go  model  which  does  not  take  this  into 
consideration  is  highly  suspect  in  performance  results. 
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(7)  No  attempt  is  made  to  optimize  the  filter  model  design  used  throughout 
this  simulation.  An  optimized  filter  design  may  (or  may  not)  provide  additional 
performance  benefits  which  need  to  be  understood  and  characterized. 

5.3  Overall  Summary 

This  thesis  investigates  using  an  extended  Kalman  filter  to  predict  a  kinetic 
kill  zone  for  aircraft  seif  defense  against  homing  missiles.  The  problem  is  studied 
using  FORTRAN  simulations  conducted  on  Sun  SPARC  II  Workstations,  and  includs  a 
software  simulation  package  called  MSOFE  [3.  The  simulation  model  is  limited  to 
constant  altitude  planar  analysis  and  assums  benign  easily  predictable  aircraft 
maneuvers,  a  kinetically  limited  missile,  and  a  bO  Hz  laser  radar  operating  in  a 
benign  turbulence  environment.  The  overall  simulation  goal  is  to  predict  a  volume 
cross-section  relative  to  the  aircraft's  future  position,  which  the  missile  would 
have  to  fly  through  at  a  specific  time.  Assumptions  and  parameters  of  the 
developed  model  are  systematically  studied  to  determine  that  parameter's  impact  on 
the  predicted  volume  cross-section.  It  is  found  that  the  filter  pseudo-noise  used  to 
drive  the  filter's  estimation  of  the  missile's  proportional  navigation  constant  has 
the  largest  impact  on  the  estimated  volume  cross-section.  Other  parameter  impacts 
are  generally  small  enough  in  comparison,  so  as  to  be  unaccessible  without  Monte- 
Carlo  simulation. 
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Appendix  A 


This  appendix  contains  simulation  plots  for  the  missile  approach  angle  run 
series.  The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to  the  north. 
All  measurement  noises  and  dynamic  noise  strengths  are  set  at  baseline  values. 
Eight  missile  shots  are  taken  at  n/4  intervals  beginning  with  0.0  radians  and 
ending  with  7x/4  radians.  Each  run  has  a  trajectory  plot,  estimated  range  and 
range  error  plots,  estimated  LOS  and  LOS  error  plots,  and  a  series  of  one  second 
filter  divergence  plots  used  to  estimate  the  volume  cross-section.  The  number  of 
volume  cross-section  estimations  made  for  each  approach  angle  vary.  The  3*/2 
approach  angle  has  two  LOS  plot  sets.  The  second  LOS  plot  is  taken  to  the  maximum 
possible  simulation  time  to  illustrate  the  zero  range  induced  filter  divergence. 
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Appendix  B 


This  appendix  contains  simulation  plots  for  the  ladar  range  noise  run  set. 

The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to  the  north.  All 
measurement  noises  (except  range)  and  dynamic  noise  strengths  are  set  at  baseline 
values.  The  missile  is  launched  from  the  n/4  approach  angle.  Each  run  has  a 
trajectory  plot,  estimated  range  errors,  range-rate  errors,  LOS  errors  and  LOS-rate 
errors.  Estimated  range  and  LOS  plots  are  found  in  the  n/4  missile  approach  angle 
run  plots  in  Appendix  A.  The  baseline  (see  Figure  4-4)  range  measurement  error  is 
0.5  meters,  with  additional  runs  made  for  range  errors  of  1.0,  2.0  and  4.0  meters. 

i 

Each  plot  set  contains  filter  divergence  plots  of  range  error  and  angle  error  for  a 
three  second  divergence  from  three  to  six  seconds.  The  divergence  plots,  which 
define  the  volume  cross-section,  are  used  for  performance  comparison  versus  the 
baseline.  In  Figures  B-2  thru  B-19  both  the  filter  and  truth  model  measurement 
noises  are  varied  from  the  baseline.  In  Figures  B-20  thru  B-37  only  the  filter 
measurement  noise  is  varied  above  baseline  (i.e.,  the  filter  noise  is  set  "high").  In 
Figures  B-38  thru  B-55  only  the  true  measurement  noise  is  varied  above  baseline 
(i.e.,  the  filter  noise  is  set  "low"). 
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Appendix  C 


This  appendix  contains  simulation  plots  for  the  ladar  range-rate  noise  run 
set.  The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to  the  north.  All 
measurement  noises  (except  range-rate)  and  dynamic  noise  strengths  are  set  at 
baseline  values.  The  missile  is  launched  from  the  n/4  approach  angle.  Each  run 
has  a  trajectory  plot,  estimated  range  errors,  range-rate  errors,  LOS  errors  and 
LOS-rate  errors.  Estimated  range  and  LOS  plots  are  found  in  the  n/4  missile 
approach  angle  run  plots  in  Appendix  A.  The  baseline  (see  Figure  4-4)  range-rate 
measurement  error  is  1 .0  m/sec,  with  additional  runs  made  for  range-rate  errors  of 
2.0,  4.0  and  8.0  m/sec.  Each  plot  set  contains  filter  divergence  plots  of  range  error 
and  angle  error  for  a  three  second  divergence  from  three  to  six  seconds.  The 
divergence  plots,  which  define  the  volume  cross-section,  are  used  for  performance 
comparison  versus  the  baseline.  In  Figures  C-2  thru  C-19  both  the  filter  and  truth 
model  measurement  noises  are  varied  from  the  baseline.  In  Figures  C-20  thru  C-37 
only  the  filter  measurement  noise  is  varied  above  baseline  (i.e.,  the  filter  noise  is 
set  "high").  In  Figures  038  thru  055  only  the  true  measurement  noise  is  varied 
above  baseline  (i.e.,  the  filter  noise  is  set  "low"). 
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Appendix  D 


This  appendix  contains  simulation  plots  for  the  ladar  angle  noise  run  set. 

The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to  the  north.  All 
measurement  noises  (except  angle)  and  dynamic  noise  strengths  are  set  at  baseline 
values.  The  missile  is  launched  from  the  x/4  approach  angle.  Each  run  has  a 
trajectory  plot,  estimated  range  errors,  range-rate  errors,  LOS  errors  and  LOS-rate 
errors.  Estimated  range  and  LOS  plots  are  found  in  the  n/4  missile  approach  angle 
run  plots  in  Appendix  A.  The  baseline  (see  Figure  4-4)  angle  measurement  error  is 
40  tirad,  with  additional  runs  made  for  angle  errors  of  80,  160  and  250  prad.  Each 
plot  set  contains  filter  divergence  plots  of  range  error  and  angle  error  for  a  three 
second  divergence  from  three  to  six  seconds.  The  divergence  plots,  which  define 
the  volume  cross-section,  are  used  for  performance  comparison  versus  the 
baseline.  In  Figures  D-2  thru  D-17  both  the  filter  and  truth  model  measurement 
noises  are  varied  from  the  baseline.  In  Figures  D-18  thru  D-33  only  the  filter 
measurement  noise  is  varied  above  baseline  (i.e.,  the  filter  noise  is  set  "high").  In 
Figures  D-34  thru  D-5 1  only  the  true  measurement  noise  is  varied  above  baseline 
(i.e.,  the  filter  noise  is  set  "low").  There  are  no  divergence  plots  for  the  "both 
varied",  and  "filter  high"  runs  at  250  nrad  because  divergence  failed  during 
simulation  and  plots  could  not  be  generated. 
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Appendix  E 


This  appendix  contains  simulation  plots  for  the  ladar  angle-rate  noise  run 
set.  The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to  the  north.  All 
measurement  noises  (except  angle-rate)  and  dynamic  noise  strengths  are  set  at 
baseline  values.  The  missile  is  launched  from  the  n/4  approach  angle.  Each  run 
has  a  trajectory  plot,  estimated  range  errors,  range-rate  errors,  LOS  errors  and 
LOS-rate  errors.  Estimated  range  and  LOS  plots  are  found  in  the  n/4  missile 
approach  angle  run  plots  in  Appendix  A.  The  baseline  (see  Figure  4-4)  angle-rate 
measurement  error  is  2.0  mrad/sec,  with  additional  runs  made  for  angle-rate 
errors  of  4.0  and  8.0  mrad/sec.  Each  plot  set  contains  filter  divergence  plots  of 
range  error  and  angle  error  for  a  three  second  divergence  from  three  to  six 
seconds.  The  divergence  plots,  which  define  the  volume  cross-section,  are  used  for 
performance  comparison  versus  the  baseline.  In  Figures  E-2  thru  E-13  both  the 
filter  and  truth  model  measurement  noises  are  varied  from  the  baseline.  In 
Figures  E-14  thru  E-25  only  the  filter  measurement  noise  is  varied  above  baseline 
(i.e„  the  filter  noise  is  set  "high").  In  Figures  E-26  thru  E-37  only  the  true 
measurement  noise  is  varied  above  baseline  (i.e.,  the  filter  noise  is  set  "low"). 
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Appendix  F 


This  appendix  contains  simulation  plots  for  the  ladar  60  Hz  and  30  Hz  pulse- 
rate  noise  run  sets.  The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn 
to  the  north.  Ail  dynamic  noise  strengths  are  set  at  baseline  values.  All  ladar 
measurements  are  varied  together  as  presented  in  Table  3-2.  The  missile  is 
launched  from  the  ji/4  approach  angle.  Each  run  has  a  trajectory  plot,  estimated 
range  errors,  range-rate  errors,  LOS  errors  and  LOS-rate  errors.  Estimated  range 
and  LOS  plots  are  found  in  the  n/4  missile  approach  angle  run  plots  in  Appendix  A. 
Each  plot  set  contains  filter  divergence  plots  of  range  error  and  angle  error  for  a 
three  second  divergence  from  three  to  six  seconds.  The  divergence  plots,  which 
define  the  volume  cross-section,  are  used  for  performance  comparison  versus  the 
baseline.  In  Figures  F-2  thru  F-25  a  60  Hz  pulse-rate  is  assumed.  In  Figures  F-26 
thru  F-49  a  30  Hz  pulse-rate  is  assumed. 
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Appendix  G 


This  appendix  contains  simulation  plots  for  the  aircraft  INS  measurement 
noise  run  sets.  The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to 
the  north.  All  dynamic  noise  strengths  and  ladar  measurement  noises  are  set  at 
baseline  values.  INS  acceleration  measurements  are  turned  off.  INS  velocity 
measurements  are  varied  as  presented  in  Table  3-4.  The  missile  is  launched  from 
the  ji/4  approach  angle.  Each  run  has  a  trajectory  plot,  estimated  range  errors, 
range-rate  errors,  LOS  errors  and  LOS- rate  errors.  Estimated  range  and  LOS  plots 
are  found  in  the  n/4  missile  approach  angle  run  plots  in  Appendix  A.  Each  plot  set 
contains  filter  divergence  plots  of  range  error  and  angle  error  for  a  three  second 
divergence  from  three  to  six  seconds.  The  divergence  plots,  which  define  the 
volume  cross-section,  are  used  for  performance  comparison  versus  the  baseline 
(Figure  4-4). 
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Appendix  H 


This  appendix  contains  simulation  plots  for  the  pilot  induced  error  run  sets. 
The  aircraft  is  initially  heading  east  and  beginning  a  3  g  turn  to  the  north.  All 
dynamic  noise  strengths  and  measurement  noises  (except  pilot  error)  are  set  at 
baseline  values.  INS  acceleration  measurements  are  turned  off.  Noise  diffussion 
strength  representing  pilot  induced  error  is  varied  from  the  (O.lg)Vsec  baseline 
value  to  (0.5g)2/sec  in  (0.1g)2/sec  increments.  The  missile  is  launched  from  the 
ji/4  approach  angle.  Each  run  has  a  trajectory  plot,  estimated  range  errors,  range- 
rate  errors,  LOS  errors  and  LOS-rate  errors.  Estimated  range  and  LOS  plots  are 
found  in  the  x/4  missile  approach  angle  run  plots  in  Appendix  A.  Each  plot  set 
contains  filter  divergence  plots  of  range  error  and  angle  error  for  a  three  second 
divergence  from  three  to  six  seconds.  The  divergence  plots,  which  define  the 
volume  cross-section,  are  used  for  performance  comparison  versus  the  baseline 
(Figure  4-4). 


H-l 


range-rate  in  meters/second 


angle-rate  in  rad/second  angle  in  rad 


range  in  meters 


FILTER  RANGE  ESTIMATION  ERROR  Pilot  Error  =  .2g 


time  in  seconds 


time  in  seconds 


H-5 


range-rate  in  meters/second 


</> 

k. 

<D 

0) 

6 


c 

0) 

o> 

c 


CO 


0,6 
0.4 
0.2 
0 

-0.2 
-0.4 

0  1  2  3  4  5  6  7 

time  in  seconds 


rILI  tH  RANGE  ERRORS  (Pilot  Error  =  ,3g) 

i - , - , - 1 - i - r 


....  true  range  error;  filter  estimated  std-dev 


w 


J - 1 _ I _ I _ I _ L 


FILTER  RANGE-RATE  ERRORS  (Pilot  Error  =  .3g) 


....  true  range-rate  error; 


filter  estimated  std-dev 


3  4 

time  in  seconds 


H-6 


angle-rate  in  rad/second  angle  jn  raCj 


FILTER  LOS  ERRORS  (Pilot  Error  =  .3g) 

— i - i - 1 - , - r 


o 

- 1 - 

r 

•  rJ' 

....  true  angle  error; _ filter  estimated  std-dev 


_L 

2 


_ i _ i 

3  4 

time  in  seconds 


j. 


i 


5  6  7 


H-7 


angle  in  radian  -je  in  meters 


FILTER  RANGE  ESTIMATION  ERROR  Pilot  Error  =  .3g 


1 00  r 


n - r 

....  true  range  error; 


.  filter  estimated  sigma 


Or - 


■1001 


time  in  seconds 


FILTER  LOS  ESTIMATION  ERROR  Pilot  Error  =  .3g 


H-8 


true  range  error;  filter  estimated  std-dev 


2 


1 _ | 

3  4 

time  in  seconds 


FILTER  RANGE-RATE  ERRORS  (Pilo 


true  range-rate  error; 


filter  estimate 


- 1 _ I _ I _ 

3  4 

time  in  seconds 


2 


angle-rate  in  rad/second  angle  in  rad 


H-10 


angle  in  radian  range  in  meters 


range-rate  in  meters/second 


0.4 


£  0.2 
0) 

E 

c  0 

0) 

D) 

I -0.2 


M 


FILTER  RANGE  ERRORS  (Pilot  Error  =  .5g) 

- 1 - 1 - 1 - 1  r 


....  true  range  error;,  filter  estimated  std-dev 


4 1 _ i _ i _ i _ i - 1 - — i - 1 

0  1  2  3  4  5  6  7 

time  in  seconds 

FILTER  RANGE-RATE  ERRORS  (Pilot  Error  =  ,5g) 

8i - 1 - 1 - 1 — -  i  i  r 


....  true  range-rate  error; _ filter  estimated  std-dev 


2 


0 


J _ I _ I - 1 - L 

2  3  4  5  6  7 

time  in  seconds 


H-12 


angle-rate  In  rad/second  angle  in  rad 


FILTER  LOS  ERRORS  (Pilot  Error  =  .5g) 


2  3  4  5 

time  in  seconds 

FILTER  LOS-RATE  ERROR  (Pilot  Error  =  .5g) 


c  i  ....  true  angle-rate  error; _ filter  estimated  std-dev 

; 

\  - 
*«,.  _ _ -  —  -  --  --  --  ■»*«••••  “xfr-“T  *  ”  -  -  - 


3  4 

time  in  seconds 


H-13 


range  in  meters 


100 

0 


FILTER  RANGE  ESTIMATION  ERROR  Pilot  Error  =  .5g 

i - 1 - 1 - 1 - r~ 


....  true  range  error; _ filter  estimated  sigma 


-1 _ I _ I _ L 

2  3  4  5 

time  in  seconds 


x  10- 


FILTER  LOS  ESTIMATION  ERROR  Pilot  Error  =  .5g 


c 

eg 

■o 

CO 


C  0 
0) 


O) 

c 

CO 


-5. 


true  error; _ filter  estimated  sigma 


0 


time  in  seconds 


H-14 


Appendix  I 


This  appendix  contains  simulation  plots  for  the  aircraft  turn-rate  and  n 
parameter  run  sets.  The  aircraft  is  initially  heading  east  and  beginning  a  turn  to 
the  north.  The  turn  rate  varies  from  1.0  g  to  5.0  g  in  1.0  g  increments.  The  n- 
parameter  noise  set  at  (0.25)/sec  and  4.0/sec  for  each  turn  (except  baseline  which 
also  includes  an  n  parameter  noise  strength  of  (1.0)/ sec).  Baseline  conditions  are  a 
3g  turn  with  n  parameter  noise  strength  set  at  (4.0)/sec  All  other  dynamic  noise 
strengths  and  measurement  noises  are  set  at  baseline  values.  INS  acceleration 
measurements  are  turned  off.  The  missile  is  launched  from  the  n/4  approach 
angle.  Each  turn  set  has  a  trajectory  plot,  estimated  range  errors,  range-rate 
errors,  LOS  errors  and  LOS-rate  errors.  Each  turn  set  contains  filter  divergence 
plots  of  range  error  and  angle  error  for  a  three  second  divergence  from  three  to 
six  seconds  for  all  values  of  n  parameter  noise.  The  divergence  plots,  which  define 
the  volume  cross-section,  are  used  for  performance  comparison  versus  the 
baseline  (Figure  4-4). 
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